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ABSTRACT 

Making use of HI 21 cm line measurements from the ALFALFA survey (a. 40) and photometry from the 
Sloan Digital Sky Survey (SDSS) and GALEX, we investigate the global scaling relations and fundamental 
planes linking stars and gas for a sample of 9417 common galaxies: the a.40-SDSS-GALEX sample. In 
addition to their HI properties derived from the ALFALFA dataset, stellar masses (Af,) and star formation rates 
(SFRs) are derived from fitting the UV-optical spectral energy distributions. 96% of the a.40-SDSS-GALEX 
galaxies belong to the blue cloud, with the average gas fraction Jhi = Mhi/M^, ~ 1.5. A transition in 
SF properties is found whereby below Af, ~ 10^ '^ A/0, the slope of the star forming sequence changes, the 
dispersion in the specific star formation rate (SSFR) distribution increases and the star formation efficiency 
(SEE) mildly increases with Af,. The evolutionary track in the SSFR-Af, diagram, as well as that in the color 
magnitude diagram are linked to the HI content; below this transition mass, the star formation is regulated 
strongly by the HI. Comparison of HI- and optically-selected samples over the same restricted volume shows 
that the Hl-selected population is less evolved and has overall higher SFR and SSFR at a given stellar mass, 
but lower SEE and extinction, suggesting either that a bottleneck exists in the HI to H2 conversion, or that the 
process of SF in the very Hl-dominated galaxies obeys an unusual, low efficiency star formation law. A trend 
is found that, for a given stellar mass, high gas fraction galaxies reside preferentially in dark matter halos with 
high spin parameters. Because it represents a full census of Hl-bearing galaxies at z ^ 0, the scaling relations 
and fundamental planes derived for the ALFALFA population can be used to assess the HI detection rate by 
future blind HI surveys and intensity mapping experiments at higher redshift. 

Subject headings: galaxies: evolution - galaxies: fundamental parameters - galaxies: ISM - galaxies: star 
formation - radio lines: galaxies - surveys 



1. INTRODUCTION 

In the last decade, the galaxy catalogs contributed by legacy 
programs like the Sloan Digital Sky Survey (SDSS) and 
the Galaxy Evolution Explorer (GALEX) satellite extragalac- 
tic surveys have enabled us to quantify properties associ- 
ated with the stellar populations of galaxies in the local uni- 
verse. Through their statistically-based insight into the stellar 
component and the interrelationship of the physical param- 
eters, these surveys have provided quantitative clues of im- 
portance to our understanding of the formation and evolution 
of galaxies (e.g. Brinchmann et al. 2004; Salim et al. 2007). 
The bi modal distribution evident in the color-mag nitude dia- 
gram (iBaldrv et al.ll2004l: ISchiminovichetani2007 ) suggests 
a likely evolution scenario whereby galaxies in the blue cloud 
form stars vigorously and grow through mergers and later, af- 
ter depleting their gas reservoirs, then migrate to the red se- 
quence. This picture is also supported by tracers of the star 
formation history (SFH): e.g., the specific star formation rate 
(SSFR), the star formation rate (SFR) per unit stellar mass, 
is see n to vary with the total stellar mass (" Brinchmann et al.l 
120041) . The star-forming sequence (high SSFR) is associated 
with actively star forming blue cloud galaxies. Indeed, the 
stellar mass appears to be the crucial quantity governing the 
star formation (SF) along this sequence. In the absence of 
mergers or other events that trigger a starburst, blue galaxies 
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on the sequence evolve towards higher stellar mass and lower 
SSFR, and eventually become red and dead. 

In statistical terms, all surveys are biased by the proper- 
ties that define them. For example, optical surveys are biased 
in terms of optical flux and possibly surface brightness. The 
SDSS legacy galaxy redshift sample has an apparent r-band 
Petrosian magnitude limit of 17.77, as well as a surface bright- 
ness limit of 23.0 m ag arcsec^^ at the half light radius in r 
(IStrauss et al.ll2002h . In contrast, blind HI surveys are unbi- 
ased by optical characteristics but have their own limitations 
in terms of HI li ne emission sensitivity, usually as a fun ction 
of HI line width (iMartin et al.ll2mOtlHavnes""etai]|20TT) . Be- 
cause Mhi /Lopt increases with decreasing Lopt, Hl-selected 
samples are more inclusive of star-forming galaxies than op- 
tical samples of similar depth. Indeed, since almost all star- 
forming galaxies contain neutral gas, an Hl-selected sample 
can approa ch a full cens us of star-forming galaxies. For ex- 
ample. Wes t et al] (120101) demonstrated the Hl-selection iden- 
tifies galaxies with lower surface brightness, smaller absolute 
magnitudes, bluer colors and smaller stellar masses than those 
used in typical SDSS studies. However, the limitation with an 
Hl-selected sample is that it will miss the early type galaxies, 
which contain very little neutral gas (iGarcia-Appadoo et al.l 
120091) . Furthermore, analysis of the spatial correlation func- 
tion ^(r) shows that the Hl-selected ga laxies repres ent the 
least clustered population on small scales (iMartin et al.l20 12). 
a fact important for interpreting the results of future HI inten- 
sity mapping experiments. 



2 



S. Huang et al. 



Exploiting the mapping capability of the Arecibo L-band 
feed array (ALFA) and the sensitivity of the Arecibo 305 m 
ante nna, the Are cibo Legacy Fast ALFA (ALFALFA) sur- 
vey (iGiovanelli e t al. 2005a b) is an ongoing blind HI sur- 
vey, aimed at mapping ^7000 square degrees of high galac- 
tic latitude sky between 0° and +36° in declination. When 
complete, the survey will detect more than 30,000 galax- 
ies out to redshift of 0.06 wit h a median recessi onal ve- 
locity, cz, of - 8200 kms^i (iHavnes et al.1 1201 Ih. Com- 
pared to the HI Parkes All Sky Survey (HIPASsT Barnes et all 
IJOOI; Meyer etal. 2004; Wong et al. 2006), ALFALFA is 8 
times more sensitive, its spectral coverage extends over 1.6 
times the bandwidth of HIPASS, and the angular resolution 
of ALFA is 4 times better than that of the Parkes multifeed 
receiver The combination of sensitivity and spectral band- 
width enables ALFALFA to detect thousands of massive HI 
disks with Mhi > 10^*^ M©. In addition, the Arecibo spec- 
tral backend yields a finer velocity resolution than charac- 
terized HIPASS, making it possible to detect sources with 
H I line widths as narr ow as ^ 15 kms^^. As discussed 
in IHavnes et al.l (1201 Ih and earlier papers in the ALFALFA 
series, the median centroiding accuracy of the HI sources is 
^ 20" allowing the identification of most probable optical 
counterparts in 97% of cases. ALFALFA provides the first 
full census of Hl-bearing objects over a cosmologically sig- 
nificant volume of the local universe dMartin et al.ll2010l) so 
that its population includes the very rare objects missing from 
the smaller volumes sampled by previous blind HI surveys. It 
enables the study of the characteristics of Hl-selected galax- 
ies in comparison with the galaxy populations included in the 
SDSS and GALEX surveys. 

The 201 1 ALFALF A catalog, 'a AO', co vers - 40% of the 
final targeted sky area (IHavnes et al.l201 lb . giving HI masses, 
systemic velocit ies a nd HI line widths for '-^ 16000 high qual- 
ity detections (fjO). In addition to the HI line parameters, 
the a. 40 catalog also includes an assignment of the most 
probable optical cou nterpart (OC) to each HI line detection. 
IHavnes et al.l ( 1201 ll) discuss the process of assigning OCs to 
the HI sources, and where the footprints of the surveys over- 
lap, the HI sour ces are cross-referenc ed to the SDSS Data Re- 
lease 7 (DR7: Abazaiian et al."2009^, permitting the deriva- 
tion of properties associated with the stellar components of 
th e Hl-bearing galaxie s detected by ALFALFA. As shown 
bv IHavnes et al.l (1201 ll) . the a. 40 population is highly biased 
against red-sequence objects. 

In this paper, we investigate further the nature of the stellar 
counterparts of the a. 40 HI line detections, adding to the op- 
tical SDSS data photometric measures from the GALEX cat- 
alog (i j2.2l i. By combining the measurements of the gaseous 
and stellar components, we characterize the a. 40 galaxy pop- 
ulation and study in f|3]global trends within it. SED-fitting to 
the seven bands from the UV to optical is applied to derive the 
principal stellar properties, including the stellar masses i^A^ 
and SF ( i)3.2t . To understand better the characteristics of the 
gas-bearing galaxies and the potential bias associated with HI- 
selection, we define a volume-limited sub-sample extracted 
from the a. 40 catalog with a similar one extracted from the 
SDSS-DR7 ( ^4Tii an d compare th e two in ^421 i n term s of 
survey depth^gXTJ, extinction ( ^4T2] i, color ( ^423[ . SF 
behavior ( i)4.2.4| i. etc. The empirical distribution of the halo 
spin parameter is derived in f|5] which suggests that the HI- 
selected galaxies favor high spin parameter halos. Our con- 
clusions are summarized in f|6] 

Throughout this paper, we adopt a reduced Hubble constant 



h = iJo/(100 kms-i Mpc-i) = 0.7. A iChabrieii (l2003l) 
IMF is adopted. 

2. SAMPLE AND DATA 
2.1. ALFALFA parent sample 

The HI parent sample used here is dra wn from the a. 40 
catalog presented in Havnes et al.l (1201 ll) . The a. 40 cata- 
log covers two regions in the Spring sky (i.e., the Virgo di- 
rection, 7^30" < RA < 16'^30'", 4° < Dec < 16° and 
24° < Dec < 32°) and two in the the Fall sky (i.e., the anti- 
Virgo direction, 22"^ < RA < 3^, 14 ° < Dec < 1 6° an d 
24° < Dec < 28°). As discussed in IHavnes et all (120111) . 
ALFALFA HI line detections are categorized according to 
their signal-to-noise ratio (S/N) and corresponding reliabil- 
ity. Code 1 sources have S/N > 6.5 and are highly reliable. 
Another set of entries, designated as "code 2" sources, or "pri- 
ors", have lower S/N but coincide with a likely OC at the 
same redshift . Most of these sources are likely to be real 
(IHavnes et al.ll201 1). The catalog includes 15041 extragalac- 
tic HI sources, 11941 with code 1 and 3100 with code 2. In 
this paper, we consider both the code 1 and 2 a. 40 detections 
and adopt the HI measures, distanc es and HI m asses presented 
in the a. 40 presented by Havnes et aLj (1201 ll) . It is important 
to note that ~' 70% of the ALFALFA sources are new HI de- 
tections proving that previous targeted surveys based on opti- 
cal selection (magnitude, size and morphology), most notably 
the exte nsive collection co ntained in the Cornell digital HI 
archive ("Springob et al.l l20b5l) . missed large segments of the 
local gas-bearing population. 

In particular, one of the most surprising results of AL- 
FALFA to date is the richness of the high HI mass galaxy 
population (Martin et al. 2010). With its combination of sen- 
sitivity and depth, ALFALFA reveals that there still exists at 
z ^ a population of massive galaxies which retain mas- 
sive HI (Mhi > 10^° Mq) disks. Some, in fact, contain a 
dominant fraction of their baryons in HI gas. As the most HI 
massive local galaxies, they ar e the z ^ analogs o f the mas- 
sive disks detected at z ^ 0.2 ( Catinella et alj2008l) and those 
that will dominate the deep surveys being planned for even 
higher z with the EVLA, APERTIF, ASKAP and MeerKAT 
and eventually the Square Kilometre Array (SKA). By design, 
ALFALFA provides a census of HI bearing galaxies within 
a cosmologically significant volume over a wider dynamic 
range of HI masses than previous studies. It thus serves as 
the reference z ^ Hl-selected population. 

2.2. Optical and UV counterparts of ALFALFA HI sources 

In addition to the HI measurements. lHavnes et all (1201 Ih at- 
tempt to identity the most probable optical counterpart (OC) 
of each HI line source in the a. 40 catalog. Ancillary informa- 
tion such as redshift coincidence, angular size, color and mor- 
ph ology is used in ma king the OC assignment. As discussed 
bv lHavnes et al.l(l201 1 ). the process is not perfect, particularly 
for low signal-to-noise sources for which the centroiding ac- 
curacy can exceed 30" and in regions of source confusion. 
Nonetheless, the vast majority of OC assignments are prob- 
ably valid and thus permit the comparison of the stellar and 
gaseous components of the ALFALFA population. 

2.2.1. SDSS photometry 

Towards this aim, 'Haynes et al.| (1201 ll) also provide a 
cross reference of the assigned OCs with the SDSS-DR7 
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dAbazaiian et al.ll2009h where the two surveys have overlap- 
ping sky footprints. The northern Fall region is not covered 
by the SDSS legacy imaging survey DR7. Of the 15041 HI 
so urce included in the a. 40 catalog, 201 have no OC assigned 
by iHaynes et all (l20i 1 ). 2310 He outside the DR7 footprint, 
and 60 appear to be in the region of the SDSS imaging survey 
but cannot be associated with an object in the SDSS photo- 
metric database. In most of the latter cases, the OC is evident 
in the SDSS images but is projected close to a bright fore- 
ground star or contaminated by its glare. 

As discussed by numerous authors e.g.. lWest et al.l (|2010() . 
gas rich nearby galaxies are often blue and patchy with the 
result that their overall optical emission is shredded among 
several photometric objects by the SDSS pipeline measure- 
ments. The ALFALFA -SDSS cross-reference given in Table 3 
of Ha ynes et al.l (1201 ll) includes a photometric code that iden- 
tifies objects with suspicious SDSS photometry, and we ex- 
clude such objects (293/12470) from the SED-fitting. For the 
others, we require that model magnitudes are available in all 
five SDSS bands (u, g, r, i, z) and retrieve the SDSS photo- 
metric parameters from the standard SDSS-DR7 database. 

The current a. 40 catalog used here relies on the opti- 
cal identifications and photometry derived fro m the SDSS- 
pR7 d at abase as cross - refere nced in Table 3 of IHaynes et alJ 
(120111) . iBlanton et al.l ( 1201 lb have presented a new method 
for background subtracting the SDSS imaging which they ap- 
ply to the SDSS-III DR8 images. As ALFALFA progresses, 
we will migrate to use of the improved SDSS pipeline. Es- 
pecially for purposes of comparison with results obtained by 
other authors who have likewise relied on DR7, we retain the 
use of DR7 in the present work. 

2.2.2. GALEX photometry 

Using a similar approach, we have conducted a separate 
cross-match of the ALFALFA OCs to the GALEX UV photo- 
metric catalog. The imaging mode of the GALEX instrument 
surveys the sky simultaneously in two broad bands, one in 
the FUV (effective wavelength of 15 16 A) and a second in the 
NUV (effective wavelength of 2267 A). The GALE X field of 
view is ~ 1.2° in diameter (iMorrissey et al.ll2007l) . although 
image quality deteriorates in the outer annulus beyond a ra- 
dius of 0.55°. We use the GALEX GR6 data release, with its 
improvements to flat fielding, adjustment to the photometric 
zero-point, etc. Given the poorer image resolution (^ 4.5" 
FWHM), compared to that of the SDSS, as well as the lower 
UV source density, visual inspection shows that the GALEX 
DR6 pipeline measurements corresponding to the ALFALFA 
galaxies suffer less from shredding issues. We have also vali- 
dated that the GALEX pipeline photometry is in close agree- 
ment with our own photometric reprocessing of the GALEX 
images for a sample of dwarf galaxies detected by ALFALFA 
(lHuangetal.il2012D . 

The GALEX mission includes several survey modes that 
differ in their exposure time per tile: The All-sky Imaging 
Survey(AIS) is the shallowest (^100 sec), while the Medium 
Imaging Survey (MIS, ^1500 sec) is designed to maximize 
the coverage of the sky that is included in the SDSS. The latter 
also includes the Nearby Galaxy Survey (NGS), a selected set 
of targeted fields of similar depth as the standard MIS. A third 
Deep Imaging Survey (DIS) is much deeper but covers only a 
small solid angle. In adopting the GALEX counterparts, we 
give preference to sources extracted from the MIS, NGS or 
DIS, but make use of the AIS for those objects not included 



in the deeper surveys. 

To cross-match the ALFALFA OCs to the GALEX cata- 
log, we first search the position of the ALFALFA OCs for all 
GALEX neighbors within 36". Objects close to the GALEX 
field edge, i.e. with a distance from the field center > 0.55°, 
are dropped to avoid duplications with objects in overlapping 
tiles and known GALEX imaging artifact effects. To permit 
homogeneous SED-fitting to all 7 bands (FUV, NUV, ugriz), 
we require that UV images of comparable depth must exist for 
both of the UV bands; this criterion results in the adoption of 
AIS measurements where the MIS catalog is incomplete, e.g. 
because of the failure of the FUV detector. In other cases, 
matches are missed because the FUV and NUV sources may 
fail to coincide because of, e.g. astrometry error or shred- 
ding in one or both bands. Among the 14840 extragalactic 
ALFALFA sources with OCs 1828 (12.3%) have no GALEX 
counterpart returned within 36", either because they lie out- 
side of the GALEX footprint or have no detected UV emission 
in both bands. 516 (3.5%) are excluded because they lie too 
close to a GALEX field edge and 1317 (8.9%) are not matched 
because all neighbors are detected only in one band but not the 
other. The remaining 11179 OCs are matched to the nearest 
neighbor in the GALEX catalog with 7752 (52.2%) matched 
to UV sources found in the AIS and 3427 (23.1%) to ones in 
the MIS. The median separation between the coordinates of 
the OC and the cross-matched GALEX object is only 1.6". 

Although the addition of the GALEX UV photometry 
sets bet ter constr aints on the SF properties and dust extinc- 
tion (Sali m et"an 1200 5). the requirement that UV sources 
must be detected in both bands introduces an additional bias 
against non star-forming galax ies. SED-fitting to the SDSS 
bands only dHuang et alj|2012h demonstrates that 13.5% of 
the 12156 galaxies in the a.40-SDSS DR7 cross-match have 
log SSFR < 1G~^^ yr~^, and therefore belong to the quies- 
cent population. This fraction drops to 3.9% of the 9417 a. 40 
galaxies that have counterparts in both the SDSS and GALEX 
photometric catalogs. As we discuss in ^ the Hl-selection 
produces the stronger bias against the red sequence. 

2.2.3. UV-to-optical colors of ALFALFA galaxies 

Figure [T] presents the optical-UV color-color diagram de- 
rived for the 9417 galaxies which have complete entries in all 
three of the ALFALFA-a.40, the SDSS-DR7 and GALEX- 
DR6 catalogs and for which the 7-band SED-fitting pro- 
duces a valid result (see |3]), denoted as the a.40-SDSS- 
GALEX sample hereafter Contours and points depict the 
distribution for the a.40-SDSS-GALEX sample in high and 
low number density regions respectively, with typical error 
bars shown in the upper left corner (pipeline magnitude er- 
rors). Colors hereafter are all correcte d for G alactic extinc- 
tion but not internal extinction (see i ]4.2.2l i. In the opti- 
cal, the SDSS pipeline extinction-corrected values are used, 
while at UV wavelengths , we adopt the E(B — V) values 
based on the maps o f Schlegel et all (119981) . the lCardelli et alJ 
(119891) extinction law with Ry = Av/E{B - V) = 3.1, 
and A{X)IE{B - V) = 8.24 for t he FUV and 8.2 for the 
NUV, foflowing IWvder et all (l2007h . Because of the band- 
width limit of ALFALFA, only low redshift galaxies with 
small K-correction are included, and we ignore this term in 
computing colors, i.e., no K-correction is applied. 

Because it contrasts the recent SF, as indicated by the UV 
light, with the total past SF, as indicated by the optical light, 
the UV-to-optical color is a stronger diagnostic of SFH than 
colors derived from the optical bands only. This result is 
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also evident in Figure [T] Based on the distri b ution of their 
GALEX-SDSS matched catalog, ISaUm et al.1 ( l2007h define 
their blue cloud galaxies as those with (NUV—r) < 4. Using 
the same cutoff, 96% of the a.40-SDSS-GALEX galaxies He 
on the blue side, suggesting that Hl-selection induces a strong 
preference for blue star-forming galaxies, or conversely, a 
strong bias against red sequence galaxies. Bluewards of this 
division, the two colors are well correlated with a slope of 
S{u — r)/5( NUV — r) ~ 0.6, w hich is comparable to what 
was found bv lWyder et all dl^, 5{u - r) /5{NUV r) 
0.5. Although only a small population (411) of galaxies ap- 
pear redward of [NUV — r) — A, the (u — r) colors of the 
red objects increase less quickly with (NUV — r) and the 
distribution of the reddest tail is nearly flat. As discussed in 
i j4.2.2l the degeneracy of (u — r) among the red populations 
is even more pronounced in an optically-selected sample with 
more red galaxies. Therefore, the SED-fitting to the SDSS 
bands only is sufficient to constrain the SF for the Hl-sele cted 
blue galaxy population in general (e.g. lHuang et al.ll2012l) but 
gives systematic overestimates of SFRs for optically-selected 
red galaxies such as those in the Virgo cl uster known to have 
quenched SF (e.g.' Hallenbeck et al.ll2012 ). i.e., it is crucial to 
include the UV bands in the SED-fitting to infer the SF of the 
red population. As a result, we adopt the UV-optical color 
(NUV — r) rather than the optical-only (u — r) in the analysis 
of SF and gas properties below. 

3. GLOBAL PROPERTIES OF THE ALFALFA GALAXY POPULATION 

To derive the global properties of the stellar components 
of the ALFALFA galaxies, we adopt the methodology of 
ISalim et aTl (120071) . In particular, stellar masses and SFRs 
are derived from SED-fitting the seven GALEX/SDSS bands. 
Further details of the method and fitting quality as applied 
to the a. 40 sample are found in iHuang et al.l (120121) which 
focuses on the lowest HI mass dwarf population. In addi- 
tion, the Gaussian prior distribution of the effective optical 
depth in V band , ry, i s applied, with the mean predicted by 
iGiovanelli et al.l (Il997h and a standard deviation of 0.55 dex. 
Such an improvement reduces the overestimate of internal ex- 
tin ction and SF R with decr easing stellar mass, as identified 
by ISaUm etall (2007) (see ^43^ for more details), but still 
accounts for the effect of dust in disk systems. In this sec- 
tion, we discuss the results for the full a.40-SDSS-GALEX 
overlap sample (9417 galaxies). 

3.L Gas and stars 

Current understanding interprets the standard SDSS color 
magnitude diagram (CMD) in terms of an evolutionary sce- 
nario under which galaxies migrate from the blue cloud to the 
red sequence as they assemble their mass. This picture is fur- 
ther reinforced by the presence of the star-forming sequence 
in the SSFR vs. stellar mass diagram; more massive galaxies 
show lower SSFRs. Consistent with this picture, one would 
expect galaxies to grow increasingly gas poor and thus hav- 
ing lower HI fractions fni (defined throughout this work as 
Jhi = Mhi /M») as they assemble their mass. Therefore, 
blue galaxies with high gas fractions indicate disks Which are 
stable against co llapse, making their SF much less efficient 
(iWest et al.ll2009l) . 

3.1.1. HI versus stellar mass 

In the last decades, many studies have investigated how the 
HI content varies with stellar properties in galaxies, such as 



morphology, luminosity, size and SF activity (iGavaz zi et al. | 
1996; Bosellietal. 2001; Kannappan 2004; Disne y et al.1 
2008; Garcia- Appadoo et alj i2009^ ' .West et al.i .2009i. 120 lOt 
Toribio et al. 201 1). Despite the complex interplay of dynam- 
ics, SF, chemical enrichment and feedback etc., the stellar and 
HI components, as well as the dark matter halo, exhibit cor- 
relations with each other. However, many of these studies 
have relied on relatively small and/or inhomogeneous sam- 
ples limited to the very nearby universe. Although the main 
scaling relations were known, constraints on the accuracy of 
these relationships, as well as the quantification of their scat- 
ter are still not well determined. Based on an Ha narrow-band 
imaging survey of ^ 400 galaxies selected from ALFALFA, 
IGavazzi et al.l (l2012allbl) also investigate the relationships be- 
tween HI and newly-formed stars, emphasizing the study of 
environment effects. Here we focus on the nature of the pop- 
ulation detected by the ALFALFA survey. 

Figure|2]illustrates the relationships of the HI mass and Jhi 
(the vertical axes) with the stellar mass and color (horizontal 
axes). The contours and points outline the distributions of 
the galaxies in the q:.40-SDSS-GALEX sample; the blue di- 
amonds and solid lines trace the average values (log y) and 
(log a;) in bins of \ogx, with a bin size of 0.5 dex in panels 
(a), (c) and (d). The number of galaxies in each stellar mass 
or (NUV—r) color bin is given at the bottom of panels (c, d). 
Typical error bars of individual galaxies are shown in the cor- 
ners of panels (a) and (d). The Spearman's rank correlation 
coefficients of the relation, rs, are shown in the upper right 
corner of all panels. Compared to similar studies that have 
previ ously probed the glob a l scaling relation s involving Mhi 
(e.g. iBofliweU et all l2009t ICatinella et all |2010), the a.40 
sample offers a more complete statistical sampl ing of the full 
range of HI and stellar masses. As discussed in iHavnes et al.l 
(1201 ll) . alfalfa's combination of sensitivity, sky coverage 
and bandwidth yields a sample that probes a wide dynami- 
cal range in HI mass (7-11 dex with a mean of 9.56 dex), 
from the most massive giant spirals with Mhi > 10^" Mq 
to the lowest HI mass dwarfs with log Mhi < 10^'^ Mq. 
In fact, the stellar mass range that is probed is slightly wider: 
6-1 1.5 dex, with a mean of 9.43 dex. As an Hl-selected sam- 
ple, a.40-SDSS-GALEX demonstrates the ability to recover 
galaxies with small 

Figure |2j a) shows the distribution of Mhi with M^,. The 
cyan dash-dotted line traces the linear fit to th e GASS sample 
of hig h stellar mass galaxies (Af, > 10^° M^; ICatinella et akl 



\og{MHi{M^)) = 0.02 log A/, + 9.52. 

Note that those authors chose to calculate log(Af^f7) rather 
than (log Mhi) because the (log Mhi) value is depressed 
by the contribution of gas-poor galaxies in their A/, -selected 
sample; a similar effect results in their adoption of \og{fHi) 
rather than (log Jhi)- In contrast, an Hl-selected sample such 
as ours does not sample the low HI fraction massive objects 
so that, as a function of (log Af*), {\og Mhi) and {log f hi) 
adequately trace the main distribution. Moreover, as pointed 
out in Cortese et al. (201 1), the distribution of fni is closer to 
log-normal than Gaussian, and thus they also prefer (log fni) 
to \og{f hi). 

We confirm the previous findings that Mhi increases with 
Af,. However, the correlation does not appear to be a simple 
linear one, i.e. 6Mhi /SAI^, is smaller at the high mass end. 
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The linear fit to the blue diamonds in Figure |2a) is 

3.117, logM, < 9; 



(logAf 



Hit 



0.712(logA4 
0.276 (log 



7.042, logM, > 9. 



(1) 



This trend is consistent with the idea that once AGNs are 
turned on in massive galaxies, gas is lost due to AGN feed- 
back. The fact that Jhi is lower in massive SF/AGN compos- 
ites than in purely SF galaxies of the same mass may be the 
cause of a similar break in slope of the star-forming sequence 
(see i)3.2| at a slightly higher transition mass, log A/, ^ 9.5). 
Furthermore, compared to the high stellar mass GASS galax- 
ies, the ALFALFA population is overall more gas-rich for the 
same stellar mass (log Af, > 10) and traces a steeper slope in 
the Mhi vs. Af, scaling relation, i.e. there is a systematically 
larger discrepancy in the typical HI content of the ALFALFA 
and GASS populations in the largest A/* bins. Besides the 
change in slope, there appears to be an increased scatter in the 
Mhi distribution below log A/, ^ 9, a regime only p oorly 
sampled by other studies. In fact, iHuang et al.l (120121) point 
out that at the lowest HI masses, ALFALFA detects a pop- 
ulation of dwarf galaxies with low fni for their Af,; some 
of these objects are dwarf ellipticals/spheroidals (dE/dSph) 
galaxies in the Virgo cluster an d may have accreted the ir cur- 
rent gas supply only recently (iHallenbeck et al.ll2012h . The 
HI gas can be easily removed in low mass systems due to 
their shallow potential wells, so that the galaxy migrates onto 
the red sequence as its SF quenches. 

Figure |2|c) shows how the HI fraction fni depends on 
Af,. The cyan dash-dotted lin e again traces the GA SS re- 
sult for the high stellar masses jCatinella et al.ll20T^ . while 
the green (upper), red (lower) and yellow (middle) dashed 
lines trace the separate samples of Hl-normal galaxie s, ones 
in Virgo and outside-Virgo respectively from Corte se et al.l 
(j201 1) who looked for trends among galaxies in different en- 
vironments. Again, the known trend that higher Af* galaxies 
have lower fni is clearly evident, with a correlation coeffi- 
cient rs — —0.85. A/» depends more strongly on fni than 
on Mhi (fs = 0.71) partly because the same measure of 
the Mf, enters also in the computation of Jhi- Compared to 
other findings, the ALFALFA population uniformly includes 
galaxies which are more gas rich for a given A/*. Their ex- 
traordinarily high Jhi indicates little integrated past SF, while 
their blue colors may be attributed to a SFH that steadily 
rises to the present day or a truly young stellar component 
(iGarcia-Appadoo et a l. 2009). 

Both the GASS and Cortese'etal] ( 1201 Ih samples include 
galaxies that have lower Jhi and lie below the ALFALFA 
HI detection threshold. For example, the Virgo cluster is 
well known to contain a significant population of HI defi- 
cient galaxies (Davies & Lewis] 119731: iGiovanelli & HavnesI 
[1985; Solanes et al. 2002) whose HI line flux densities are too 
low for them to be detected by the short ALFALFA observa- 
tions; their detections were made using longer duration, target 
Arecibo observations. The offset of the ALFALFA population 
from the other samples is therefore as expected. However, it is 
interesting t o note that the scalin g of fni with log Af, derived 
here and by ICortese et al.l (1201 1) for the Hl-normal galaxies, 
while they do not coincide in amplitude, do show comparable 
slopes at intermediate masses, and perhaps the same is true 
for all samples at log A/* < 9.7. A "fast", shallow survey like 
ALFALFA derives the same trend as one which relies largely 
on much deeper, pointed observations. The flattening off of 
Jhi at log A/* < 9 is traced only by the ALFALFA dwarfs 



( IHuang et al.ll20T2 t). 

It is important to note that the nearby, low mass galaxies are 
the ones most susceptible to shredding by the SDSS pipeline 
so that, statistically, their stellar masses are more likely to be 
underestimated, resulting in an extreme tail of galaxies with 
unrealistically high fni- By exclusion of objects with sus- 
pect SDSS photometry as noted by Ha vnes et al.l (1201 ll) . the 
most egregious cases have been excluded from this analysis. 
Similar problems with the use of the SDSS pipeline measure- 
ments have much less effect on the main distribution. At the 
same time, source confusion within the ALFA beam (FWHM 
^ 3'. 5) is more likely among more distant systems so that 
the Mhi (and Jhi) of some high Mhi sources may be over- 
estimated. However, other than cases of major mergers, the 
highest Mhi galaxies are always significantly more massive 
than their small companions, so that the change in Jhi, if 
the contribution from companions is removed, would only 
be small. Overall, the trend of falling Jhi with increasing 
A/, seen in the ALFALFA galaxies is well defined. For the 
ALFALFA population overall, the median Jhi ^ 15. HI- 
selected galaxies are uniformly gas rich for their stellar mass 
following a scaling relation over the range of stellar mass 
8.0 < log A/, < 11.0. 

In Figure |2jc), the number density of points drops sharply 
on the upper edge of the main distribution: there is a real 
cutoff in the galaxy population with even higher Jhi than 
ALFALFA detects. The increased dispersion in the contours 
on the lower edge of the Jhi distribution with substantial 
numbers of outliers with lower Jhi than the main population 
confirms that, because of its Hl-selection, ALFALFA misses 
much of the gas-poor galaxy population. Longer integration 
times would obviously detect galaxies of lower Mhi and thus 
lower Jht at consta nt Af*. The GASS observing strategy 
(ICatinella et al.l2 01^ is specifically designed to probe to con- 
stant Jhi by conducting significantly longer but targeted HI 
observations. The GASS program thus characterizes the over- 
all population of galaxies selected by stellar mass at the high 
mass end. The ALFALFA survey, on the other hand, samples 
well the full dynamical range of the Hl-rich (for their stellar 
mass) population 

Figures |2b, d) explore the variation of Mhi and Jhi with 
{NUV — r). Definitions of diamonds and lines in panel 
(d) are the same as in panel (c). As noted in Figure [T] 
nearly all ALFALFA galaxies are blue, and the population 
is highly biased against the red sequence. While there is a 
wide spread in Mhi, there is little trend of Mhi with color 
(rs = 0.28). In fact, there are 128 a.40-SDSS-GALEX red 
galaxies [NUV — r > 4) with log Mhi > 10, including the 
early type galaxies with quenched SF but unusually high HI 
masses (e.g. AGC 260442), the edge on galaxies with sig- 
nificant internal extinction (e.g. UGC 6312 has a dust lane 
evident in the SDSS ima ge), and even th e "red spirals" found 
in the Galaxy Zoo ( Masters et al.l l2010). e.g. UGC 9624 and 
UGC 9283. There are 116 red galaxies (NUV - r > 4) 
with logA^g j < 9 .5; most are early type "dead" galaxies. 
ICortese et al.l (1201 ll) found that the blue cloud galaxies have 
the same Jhi regardless of environment, whereas for the red 
galaxies, Virgo members are significantly gas poorer than Hl- 
normal systems (see dashed lines in Figures|2}l). In particular, 
their fit for Hl-normal galaxies (green) agrees well with the 
main trend for the a. 40 galaxies. 

In contrast, Jhi is a strong function of color among the 
ALFALFA population (rs — —0.79 in panel d), at least 
among the blue cloud galaxies, i.e. the bluer galaxies tend 
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to have higher HI fractions. However, this trend gradually 
flattens for the very red galaxies (NUV — r > 3.5), i.e. 
the very red galaxies among the ALFALFA population have 
higher fni than would be predicted by extrapolation of the 
trend traced by the blue galaxies. Comp ared to the f nr 
vs. (NUV - r) tre nds derived by Catinella etal] (1 2010!) or 
ICortese et al.l (1201 ll) . traced by the dash-dotted and dashed 
curves respectively in panel (d), the offset of blue diamonds 
on the blue side is small, but the deviation becomes systemat- 
ically larger in the redder bins. Such a change in behavior can 
be partly explained by the fact that ALFALFA detects only a 
very small subset of these red galaxies. The presence of HI 
in this small population of otherwise "red and dead" galax- 
ies is most easily explained if their HI gas has been acquired 
only recently, as has been invoked previously to expl ain the 
HI in ellipticals (e.g., Wardle & Knapp 1986; Morgant Tet all 
2006h. and t he annular HI distributions seen in many SOs (e.g.. 



van Driel et al. 1988: Donovan et al. 2009). Deep HI synthe- 



sis imaging of the SAURON and ATLAS^^ samples of early- 
type galaxies shows that HI is commonly de tected in galax- 
ies which do not resid e in cluster cores (e.g.. lOosterloo et all 
l2010t[Serra et al.|[2012l) . A significant fraction of non-cluster 
early-types contain some cool HI gas, with the large spread 
in HI content likely due to differences in their accretion histo- 
ries. 

3.1.2. Predictors of HI gas fraction 

The tight correlation between fni and (NUV — r) can be 
used as a predictor of Mhi given meas ures of color, i.e ., 
the 'photometric gas fraction' technique (lKannappanll2004l) . 
Furthermore, a 'fundamental plane' of fnt -jNUV — r)- 
jjL-t ha s been identified by the GASS survey (ICatinella et al.l 
120101) . where the stellar mass surface density is defined as 
/X, [Mq kpc-2] = O.SAf, [Mo]/7r(r5o,2 [kpc])^ and r^o^^ is 
the radius containing 50% of the Petrosian flux in the z-band. 
Their best fit 'plane' is 

log !hi = -0.240(iV[/y - r) - 0.332 log + 2.856, 

and the scatter of such a fui predictor is reduced relative to 
the fHi-{NUV~r) correlation with the additional parameter 



/i*. Because colors essentially trace the SSFRs (see also ^3.2\ , 
similar predictors are calibrated bv lZhang et al.l (l2009h for an 
optically-selected sample as: 



log fni = 



-1.25(5 -r)- 0.54 log + 4.66; 
0.26 log SSFR - 0.77 log /i* +8.53. 



It should be noted that no correction for internal extinction is 
applied in those analyses, although it is wel l known that the 
inner disks of spirals are optically thick (e.g.. lGiovanelli et~an 
Il995h . We dis cuss t he need for an internal extinction correc- 
tion below in i|4.2.2| The fni predicted by these formulae are 
plotted on the horizontal axes in Figures [3|a-c), respectively. 
Compared to the ALFALFA measurements of the fni (ver- 
tical axes), they all predict systematically smaller Jhi- This 
reaffirms that the Hl-selected sample is biased towards the gas 
rich population. The deviation from the one-to-one dashed 
line increases wit h fni in the case of the GASS calibration 
(Figure [3^). The iZhang et al. (2009) estimators (Figure [3]5 
and c) systematically underpredict fni of the a.40-GALEX- 
SDSS galaxies by -0.3 dex. 

Exactly which scaling relation to use to predict the prop- 
erties of a population depends of course on what the scien- 
tific objective is. For example, the scaling relations for an 



optically-selected sample may be valid for a stellar mass se- 
lected population. However, the relations derived for AL- 
FALFA give better predictions for the HI detection rate for 
future SKA surveys which will likewise be Hl-selected. Pre- 
vious simulations of SKA detection rate are mostly based on 
the HIPASS results loca lly (e.g., A bdalla & Rawlings 2005]; 
lObreschkow et al.ll2009l) . But HIPASS suffered from limita- 
tions in its volume sensitivity, and in fac t, ALFALFA detects 
more HI sources at the high Mhi end dMartin et al.l 12010). 
Based on the a.40-GALEX-SDSS galaxies, linear regression 
gives: 

log Jhi = -0.25(7VC/y-r) -0.57 log M* +5.24; (2) 
= -1.05(5 -r) -0.571og/i* +5.12; (3) 
= 0.271ogS'S'Fi?- 0.64 log + 7.80. (4) 

Note that for comparison with other authors, the colors used 
here are corrected for Galactic extinction but not for inter- 
nal extinction. In Figures Od-f), the ALFALFA fni mea- 
sures are plotted against the values predicted by equations 
(|2|-(|4l), with the correlation coefficients rs = 0.88, 0.87 
and 0.87 respectively. The systematic offset is removed ac- 
cording to our best fit and the correlations are as tight as 
the Catinella et al. (2010) (rs 0.88) and the two different 
IZhang et al. (2009) (rs = 0.87 and 0.86) results. Among the 
three planes, the fHi-{NUV —r)-ii^, correlation has the least 
scatter and it is also tighter than the fHi-{NUV — r) corre- 
lation (rs = -0.79). 

Note that the fundamental plane found here is noticeably 
different from the GASS relationship (Catinella et al. 2010). 
The main trend in Figure Od) reveals a break in slope at 
log//// ^ —0.5, whereas the GASS relation is confined to 
only below this critical ////. Although the GASS sample is 
complete in the massive il/* domain, it does not probe the 
lower stellar mass, gas rich systems. Blindly applying the 
fundamental plane defined by the massive gas-poor galax- 
ies through extrapolation into the gas-rich regime results in 
serious underprediction o f ////. Because of the change in 
slope, the deviation from (ICatinella et al.ll2OT0i) is systemat- 
ically larger in the high //// galaxies. 

Because the HI population is so overwhelmingly domi- 
nated by blue cloud dwellers and since HI is presumably 
a constituent of a galaxy's disk, not its bulge or halo, it 
seems appropriate to explore scaling relations which are 
tied more heavily to the galaxy's disk. Hence, we define 
a disk stellar mass surface density //*,r.90 [M© kpc~^] = 
0.9Af* [MQ]/7r(rgo.,- [kpc])^, where rgo.r is the radius con- 
taining 90% of the Petrosian flux in the r-band. Compared to 
Ai*, /^*,r90 is based on the r-band flux with higher S/N and less 
bulge contribution. In addition, adopting colors corrected for 
internal extinction (see ^4.2.2i . we derive improved predictors 
as follows: 

log//// = -0.17(7VC/y-r)o -0.811ogAi,,,9o + 6.3l!(5) 
= -0.70(.g-r)o -0.79 log Ai,,,9o + 6.16(6) 
= 0.22 log S-S-Fi? - 0.78 log ^^,,.90 + 8.03(7) 

The values given by equation (|5]l-(|7]i are plotted on the hori- 
zontal axes in Figures[3lg-i). They show less scatter from the 
one-to-one dashed line with better rs — 0.90. We suggest that 
scaling laws which incorporate properties which reflect the 
disk nature of the HI distribution, and specifically the above 
relations, provide the most appropriate approach to predicting 
the characteristics of Hl-selected populations. 



a. 40: the Galaxy Population 
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3.1.3. Assessing the molecular gas H2 contribution 

ALFALFA is an extragalactic HI line survey and, as such, 
probes only the neutral ISM. Yet, the process of SF in most 
galaxies is more directly coupled to the molecular gas, and 
the question of which gas component - HI, H2 or total gas 
- correlates best with SF is still debated. To account for the 
full gas content, we thus need to assess the expected contribu- 
tion of molecular gas to the total gas mass in the ALFALFA 
population galaxies. 

An outgrowth of the GASS survey of high stellar mass 
galaxies, COLD GASS is a legacy survey which has mea- 
sured the CO(l-O) line of -350 randomly selected GASS 
sample galaxies (0.025 < z < 0.05) with the IRAM 30m 
telescope. COLD GASS has uncovered the existence of sharp 
thresholds in galaxy structural parameters such as con- 
centration index and {NUV — r) color, above which the de- 
tection rate of the CO line drops suddenly. These thresholds 
correspond approximately to the transition between the blue 
cloud and red sequence (Saintonge et al. 201 lb). Even though 
ICatinella et alj (12010) found some red sequence galaxies with 
a surprisingly large HI component, none of the 68 galaxies in 
the first installment of COLD GASS with [NUV - r) > 5 
are securely detected in CO. At the same time, only 1 .4% of 
the a.40-GALEX-SDSS galaxies have {NUV - r) > 5, so 
that the Hl-selected galaxies should have a high detection rate 
in CO. 

Under the assumption that molecular gas forms out of lower 
density clouds of atomic gas, one might naively expect a tight 
correlation between Mhi and ■ However, within the sub- 
sample of galaxies detected both in HI and CO by COLD 
GASS, the fraction {Mh^/Mhi) varies greatly, from 0.037 
up to 4.09; the two q uantities are only weakly correlated 
dSaintonpe et al.ll201 lah . The relative proportions of molec- 
ular and dense atomic gas in giant molecula r clouds depend 
on the cloud column density and metallicity dKrumholz et al.l 
I2OO8I) . and the clouds could ev en be primarily at omic if the 
metallicity is sufficiently low (Ostriker et al. '20101). 

Of all the parameters that Saintonge et al. (201 13) investi- 
gated, the mean molecular gas fraction (//j^ = A/^f^/^^*) 
among the COLD GASS galaxies correlates most strongly 
with their (NUV — r) color, with 

log = -0.219(7V[/F ~r)~ 0.596, 

although it is weaker than the f hi-{NUV — r) correlation, 
probably because H2 resides in the inner region where ex- 
tinction is higher, wher eas HI dominates in the outer disks 
jSaintonge et al.l 12011 al) . At the same time, they find that 
is only a weak decreasing function of Af,. As a result, 
although ALFALFA probes a lower stellar mass range than 
COLD GASS does, the fn^-iNUV - r) correlation above 
can still be roughly applied to the a.40-GALEX-SDSS galax- 
ies. Specifically, since the a.40-GALEX-SDSS galaxies have 
a mean (NUV-r) of 2.24, the results of COLD GASS predict 
a mean ///^ of 0.082 for the Hl-selected population, higher 
than the 0.066 of the COLD GASS detections. 

We note that although /^f^ is only a weak decreasing func- 
tion of Af,, fni clearly decreases with increasing Af», i.e., 
the MH2/MH1 fraction appears to decline in less massive 
galaxies (see also Blanton & Moustakas 2009). For lumi- 
nous galaxies, a substantial fraction of the gas is sometimes in 
molecular form, but the detection of CO i n low mass galaxie s 
has been shown to be very difficult (e.g., iLerov et al.ll2009h . 
Therefore, we conclude that /i/^ ^0.1 for the ALFALFA 



population and thus ignore its contribution to the total gas 
fraction, focusing instead on the well-determined atomic gas 
fraction fni- 

3.2. Star formation properties 

In addition to the stellar mass, SED-fitting also yields an 
estimate of t he cur rent SFR averaged over the last 100 Myr. 
iSalim et akl (120051) have shown the importance of including 
the GALEX UV bands, especially the FUV, to reduce the 
uncertainties in SFRs derived from SED-fitting. In addition 
to the SFR itself, several other quantities of physical inter- 
est also become available. For example, the SSFR, defined 
as SFR/M^,, compares the current SFR with that in the past 
(as measured by A/*), and thus is well correlated with the 
birthrate-, or 6— parameter, defined as SFR/ (SFR). Both 
the SSFR and the fe— parameter describe the SFH. At the same 
time, normalization of the SFR by Mhi instead of Af, yields 
the star formation efficiency (SEE), defined as SFR/Mhi- 
The SFE compares the current SFR with its potential in the 
future, the latter regulated by Mhi, the available fuel. The 
reciprocal of the SFE is the R oberts time, tji ~ Mhi / SFR 
( lRol3ertslll96llSandapelll986l) . the timescale for depletion of 
the HI gas reservoir, assuming a constant SFR at the current 
level. 

Figure |4] shows a montage illustrating how the SF related 
properties, SFR, SSFR, SFE (y axes) vary with Af*, Mhi 
and the [NUV — r) color (x axes). As before, contours and 
points trace the a.40-SDSS-GALEX population. Spearman's 
rank correlation coefficients are shown in the lower left cor- 
ners of all panels. Typical error bars of individual galaxy 
estimates are plotted in the lower right corners of panels (a, 
e, i). In the bottom row, tracing the SF E, the cyan dashed 
line sh ows the average value obtained bv lSchiminovich et al.l 
dloTol) for the high stellar mass GASS sample, while the green 
dash-dotted line marks the value corresponding to the Hubble 
timescale. 

3.2.1. The SFR and SSFR in Hl-selected Galaxies 

Among the ALFA LFA popula tion and in agreement with 
previous studies, e.g. Salim et al. (2007), SFRs generally in- 
crease but SSFRs decrease with increasing stellar mass (Fig- 
ure 2^, d). Similar trends are also evident with Mhi (P^g- 
urel^', e), albeit with larger scatter, especially in the SSFR 
vs. Mhi distribution (rs = —0.31). The trend of decreas- 
ing SSFR with increasing stellar mass suggests the "down- 
sizing" scenario of structure formation (Cowie et al. 1996), 
such that the high A/, galaxies form most of their st ars in 
the first — 3 Gyr after their formation (iBell et al.ll2003i) and 
today exhibit relatively suppressed SF. In contrast, the low 
mass systems in such a picture remain active in SF throughout 
their histories. Under a hierarchical dark matter halo assem- 
bly scenario in which the low mass structures form first and 
then merge to form massive galaxies, the "downsizing" con- 
cept suggests a late epoch of gas replenishment and regrowth 
in low mass systems. 

Although an Hl-selected sample Uke ALFALFA is biased 
against massive galaxies with low SFRs and low SSFRs 
(see also i j4.2.4b . there is a hint in Figures Ufa, d) that the 
number density of such galaxies increases in the A/» > 
10^° Mq regime in comparison to the intermediate mass 
range (10* M© < A/» < 10^° Mq). The presence of some 
points in the lower right corner of the SFR vs. A/» plot 
suggests that at least some objects with large stellar masses 
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and detectable HI but very low SFRs are included in the AL- 
FALFA population. More importantly however, there is not 
a comparably rich population of massive HI disks with low 
SFRs, i.e., the number density of galaxies in the lower right 
corner of Figure |4j3 is lower than that in the lower right corner 
of|4^. Where there is a lot of HI, there is always some SF. 

As evident in Figure SJc), the expected correlation be- 
tween SFR and {NUV — r), that galaxies bluewards of 
(NUV — r) 4 have higher S FRs than ones redwards of 
that value (e.g. Salim et al.l l2005h . is not so well defined by 
the ALFALFA population (rg = 0.31), mainly because AL- 
FALFA detects only a few very red galaxies. In particular, we 
lack sufficient dynamic range in (NUV—r) to probe the trend 
along the red sequence seen in optically-selected samples that 
galaxies bluer in (NUV — r) have higher SFRs, especially if 
colors after extinction correction are plotted. 

On the other hand, the SSFR is a strong function of the 
(NUV — r) color (Figure |4f), with the Spearman's rank cor- 
relation coefficient rs = —0.76. This is also the tightest 
among all the correlations shown in Figure |4] It is natural 
to expect that (NUV - r) is closely tied to the SSFR. Since 
the NUV luminosity largely reflects the SF and the r-band 
luminosity the stellar mass, the (NUV — r) color, as the ratio 
of the two, serves as a proxy for the SSFR. Given the fact 
that NUV better characterizes the SFR than the w-band, one 
may also expect the SSFR to correlate better with (NUV—r) 
color than with (u — r). However, because NUV suffers more 
from internal extinction and the associated corrections can be 
highly uncertain, extra scatter is introduced when the NUV 
is used with accounti ng for the impact of extinction. In fact, 
as demonstrated in ( ii4 2.2l l and Figure [HJb), a shift towards 
blue colors and an even tighter correlation between SSFR 
and (NUV — r) become apparent when extinction-corrected 
colors are used. 

3.2.2. The star formation taw in Hl-selected galaxies 

The underlying question linking gas to stars in galaxies, 
the "star formation law" (SFL), is what limits SF: the forma- 
tion of molecular gas out of HI or the efficienc y at which the 
availa ble molecular gas is converted into stars (ISchruba et al.l 
1201 IV Various forms of the SFL are studied, perhaps most 
common among them the empirical law describing how the 
SF surface density (Hsfr) i s regulated by th e gaseous sur- 
face density (e.g., T,hi+H2 inl Kennicutil 19981) . 

It should be noted that since most galaxies are unresolved 
by the ALFA 3.5' beam, ALFALFA measures only the global 
HI content. Our estimate of the SFL will thus be globally 
averaged. Numerous recent studies focusing on more de- 
tailed observations of smaller yet representative samples have 
demonstrated the regulation of SF by molecular gas. For ex- 
ample, the HI Nearby Galaxy Survey (THINGS; Walter et all 
|2008) and the HERA CO Line Extragalactic Survey (HER- 
ACLES ;|LeroiLeLalJ|2Q09) provide measurements of the sur- 
face densities of total gas, atomic and molecular gas, and SFR 
in ~kpc-sized regions within a number of nearby galaxies. 
Measurements of the azimuthally averaged gas and SFR pro- 
files show that the SFR correlates better with the molecular 
hydrogen component than with the total gas density within 
the optical disk (e.g. Bigiel et al. 2008), suggesting that the 
SFR is controlled by the amount of gas in gravitationally 
bound clouds and that H 2 is directly important for cooling. 
iKrumholz et al.l ( 1201 lah collated observations of the relation- 
ship between gas and SFR from resolved observations of 
Milky Way molecular clouds, from kpc-scale observations of 



Local Group galaxies, and from unresolved observations of 
both disk and starburst galaxies in the local universe and at 
high redshift. Those authors showed that the data are consis- 
tent with a simple, local, volumetric SFL and that the SFR 
is simply ^1% of the molecular gas mass per local free-fall 
time. Furthermore, Schrubaetal. (2011) found a tight and 
roughly linear relationship between IR (inferring 'Ssfr) and 
CO (inferred S/fg) intensity, with J^H2/^sfr ~ 1-8 Gyr. 
This relation does not show any notable break between re- 
gions that are dominated by molecular gas and those domi- 
nated by atomic gas, although there are galaxy-to-galaxy vari- 
ations in the sense that less massive galaxies exhibit larger 
ratios of SFR-to-CO, an effect which may due to depressed 
CO relative to H2 in low metallicity galaxies. Similarly, 
'Bigiel etaTI ( 1201 Ih demonstrated a roughly constant H2 con- 
sumption time (Y.H^/'^SFR ~ 2.35 Gyr). 

However, other works show that the relationship between 
SF and gas vari es syste r naticall y depending on the local en- 
vironment. Big iel et al.l (l20I0allbl) . found an evident correla- 
tion between SF and HI in the outer disks of spirals and in 
dwarf galaxies where HI is likely to dominate the ISM. Given 
the poor correlation between HI and SFR found inside star- 
forming disks, this finding strongly hints that different physics 
governs the formation of star-forming clouds, and that the 
HI colu mn is perhaps the ke y environmental factor in setting 
the SFR dBigiel et al.ll20I0ah . Furthermore, the SFL is Ukely 
to have a distinct form in the atomic-gas-dominated regime 
(e.g. E.s'FR oc ^Hi+H2\fPsd, theoretically by Ostrike ret alj 
2010, where psd is the midplane density of stars plus dark 
matter). Therefore, we may expect a steeper dependence of 
SsFJ? on Yjhi+H2 if there is a dropoff in the stellar and 
dark matter density with radius. There is no single univer- 
sal slope predicted for the SFL in the diffuse-gas-dominated 
regime. In low gas surface density or low metallicity regions 
where gas is significantly atomic, thermal and chemical pro- 
cesses become dominant in determining where stars can form, 
and the gravitational potential of the stars and dark matter 
m ay have a significant effec t. Similarly, the model developed 
bv lKrumholzetai] (l2009bl ) suggests S5F_R becomes a steep 
function of T,hi+H2 when complexes of gas become primar- 
ily atomic, for low ISM surface density. Observations also 
confirm steeper slopes for the low density outer Hl-dominated 
regions of spiral galaxies, as well as dwarf galaxies, com- 

S iared to the inner molecular-dominated regions of spirals 
Bigiel et al.ll2010al) . 

The increasing SFR with HI mass evident in Figure |4|b) 
suggests the regulation of SF by the HI gas, with a correlation 
coefficient of rs = 0.71. The red dash-dotted line in Figure 
IHb) shows the linear fit to the q;.40-SDSS-GALEX galaxies, 
with a slope of 1.19, suggesting a global, atomic, volumetric 
SFL. The close to unit slope indicates a SEE close to con- 
stant as a function of Mhi- The correlation between SFR 
and Mhi appears to be in conflict with the earlier finding that 
most galaxies show little or no correlation between T^sfr and 
T^Hi dBigiel et alj|2008l) . The high Jhi galaxies represented 
by the ALFALFA population appear to obey an unusual SFL 
which may not only depend on the H2, but also on the HI, stel- 
lar and dark matter properties. Additional observational and 
theoretical work is needed to evaluate how the SF efficiency 
of bound clouds depends on the relative amounts of cold HI 
versus molecular gas. 

3.2.3. The SFE in Hl-selected galaxies 
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Figures Sfg-i) illustrate the distribution of the SFE with 
M,, Mhi and the (NUV — r) color The timescale for atomic 
gas depletion for the majority of the ALFALFA galaxies is 
shorter than the Hubble time tn, and comparable to it for 
many of those with low stellar masses, A/» < 10^ Mq. 

For the high stellar mass GASS population, 
ISchiminovich et all (12010) found that, unlike the SSFR which 
decreases with increasing stellar mass, the SFE remains rela- 
tively constant with a value close to SFE = 10~^'^ yr or 
equivalently tji ^ 3 Gyr. This val ue is lo nger than the molec- 
ular gas depletion timescale (see ^3.2.2i . Furthermore, those 
authors also found little variation in the SFE with stellar mass 
surface density /i*, the (NUV — r) color or the concentration 
index, a result which they interpreted as an indication that 
external processes or feedback mechanisms which control 
the gas supply are important for regulating SF in massive 
galaxies. Considering that log SSFR ~ log SFE + log fni, 
an interesting implication of the weak correlation between 
the SFE and the stellar mass is that the fit to the log SSFR 
versus log A/* distribution would have a similar slope to 
that of the log /hi vs. logAf*, specifically —0.288 for 
l og M ^, < 9 and —0.724 for log Af* > 9 (see equations [T] in 
^3.1\ . However, the red dash-dotted line in Figure|4|d) shows 
the linear fit to the 'star-forming sequence' defined by the 
a.40-SDSS-GALEX galaxies: 



log SSFR = 



-0.149 log A/, 
-0.759 log A/, 



8.207, log A//, < 9.5 



2.402, log A//, > 9.5.' 



(8) 



The differences in the slopes suggest that (SFE) is a weak 
increasing function of (M^,) in the low M* range but re- 
mains relatively constant above logAf* ^ 9.5; this trend is 
also evident in the bottom row of panels(r5 = 0.35). The 
mild trend of increasing SFE with stellar mass seen in Figure 
4tg) w as not evident in the GASS study ( Schiminovich et all 
2010^ because the GASS sample includes only galaxies with 
AT* > 10^" Mq. 

Rather than a simple continuous scaling relation, the change 
of slope given in equation (O and evident in Figure SJd), 
suggests that a transition mass exists at Af* ^ lO^ '^ Mq 
in the way in which star formation scales with total AI^,. 
A similar t ransition ma s s at ^ lO^ '' M© in SSFR is 
adopted by ISalim et aLl (l2007h for their blue galaxies with 
{NUV — r) < 4. Those authors suggested that the lower 
SSFR is a consequence, at the high A/, end, of a population of 
systems which are both star-forming and have AGN, thereby 
yielding lower SSFRs than p ure SF galaxies of the same mass. 
Similarly, Kannapp an et al.l (|2009) identified a "threshold" 
stellar mass of several times 10^ Mq, below which the num- 
ber of blue sequence E/SO galaxies sharply rises. Those au- 
thors matched the threshold to the mass scale below which 
the mean HI content of low-z galaxies increases substantially 
both on the red sequence and within the blue cloud. Abrupt 
shifts in the SFE and gas richness near the "threshold mass" 
have been linked to the interplay of gas infall, supernova- 
driven winds, and changes in mass surface density. How- 
ever, it is important to note that such a threshold falls be- 
low the "transition" mass characteristic of the "green valley", 
identified in the, e.g., /hi versus A/* and /i* relations, at a 
stellar mass Af* ^ (2 — 3) x 10" ' '° M© proposed in many 
other works (e g. 'Kauffmann et all 120031; iBaldrv et al.ll2004t 
iBothwell etal] 12009; Catinell aeTaLl 120101) and suggested to 
be indicative of the SF quenching in massive galaxies as they 
migrate from the blue cloud to the red sequence. 



Perhaps surprisingly, the ALFALFA galaxies have on av- 
erage lower SFE, or equivalently, longer tji, compared to 
the optically-selected population, with a mean of SFE — 
2^0-9.95 yj. -I qj. equivalently tn — 8.9 Gyr, compared to the 
tfj^ 3 Gyr derived for the GASS galaxies. We note t hat the 
average t f( value was volume corrected in Schiminovich et al ] 
(120101) . but not in Figur e[4|g -i). However, we confirm that the 
volume correction (see i]4.1b results in only subtle changes in 
the mean tn as a function of A/* for the q:.40-SDSS-GALEX 
galaxies, and it i s still longer than tn ^ 3 Gyr. As we demon- 
strate in ij4.2.4[ the Hl-selected galaxies have, on average, 
higher SFRs at a fixed stellar mass, so that the lower SFEs 
must result from their higher HI masses rather than from less 
active states of SF. This result reaffirms the general conclusion 
that Hl-selected samples are strongly biased towards the most 
gas-rich galaxies. In agreement with the low SFEs charac- 
teristic of the ALFALFA population, Bigiel et al. (2008) have 
seen a decrease i n SFE in the Hl-dom inated THINGS galax- 
ies. Furthermore. iBigiel et alj (l2010al) found that the SFE de- 
creases with galactocentric radius among the THINGS sam- 
ple across the outer disks beyond the optical radius, where 
HI dominates the ISM, with tn well above Hubble time. In 
the THINGS dwarf galaxies, the contribution of H2 to the to- 
tal gas budget is generally sm all even in the inner disks, also 
corresponding to a low SEE (iBigiel et al.ll2010al) . All these 
results are consistent with the conclusion that SFEs are low, 
on average, in Hl-rich systems. 

However, we note that the low HI SFE may not be in con- 
flict with the usual H2 SFL. The Hl-selected high fni galax- 
ies may still follow the normal behavior of how stars form 
from H2, but rather that a bottleneck exists in the process 
by which star-forming molecular clouds assemble. The con- 
version of HI to H2 depends on environment inside a galaxy 
and the relative abundance of HI and H2 is key to setting the 
SFR (Bigiel et al. 2010b). Although the low HI SFE suggests 
the inefficiency of HI-to- H2 conversion, the H I-to-H2 ratio 
cannot be arbitrarily high. lOstriker et al.l (120101) assumed an 
equilibrium state, in which cooling balances heating and pres- 
sure balances gravity. This balance can be obtained by a suit- 
able division of the gas mass into star-forming (gravitational 
bound) and diffuse components such that their ratio is propor- 
tional to the vertical gravitational field. If too large a fraction 
of the total surface density is in diffuse gas, the pressure will 
be too high while the SFR will be too low. In this situation, 
the cooling would exceed heating, and mass would drop out of 
the diffuse gas component to produce additional star-forming 
gas. 

Close examination of the SFE vs. A/» diagram in Figure 
Ug) also reveals a considerable number of outliers with rel- 
atively high SFE at the high Af* end, falling well above the 
main distribution. In general, high SFEs have been measured 
in starburst galaxies and interacting systems which are con- 
suming their gas reservoirs on very short timescales. How- 
ever, a close inspection of 13 q;.40-SDSS-GALEX galaxies 
with log SFE > —8 shows that 9 of them are members of 
the Virgo Cluster. Ram pressure stripping results in strong 
HI deficiency and very short tj^ in these extreme outliers. In 
contrast, the outliers below the main distribution can either be 
red, massive, low SSFR galaxies against which ALFALFA is 
strongly biased, or abnormally gas rich (for their stellar mass) 
but quiescent ones. The latter include candidates of recent re- 
accretion or systems in which the HI gas is somehow inhibited 
from forming stars. We return to this point in S|5] 

Figures |4jh, i) illustrate that the SFE barely changes with 
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either Mhi or color {rs = 0.03 and 0.18 respectively) and 
the scatter in both correlations is large. 

3.2.4. Linking the gas fraction fni to SF 

As first discussed bv lRobertsI d 1 963h . it is not surprising that 
SF appears to be regulated by gas content. We have already 
argued that the log SSFR vs. log Af, diagram (Figure |4}l) 
is similar to the log fni vs. log one (Figure both 
showing similar slopes along the main trend; the distributions 
of log SSFR and log fni in given stellar mass bins also be- 
come broader at both the high- and low- mass ends. Hence, 
the star-forming sequence in the former diagram can also be 
understood as a sequence o f gas- depletion in the latter one. 

Previously, iKannappanI (120041) also linked the fni to bi- 
modalities in galaxy properties. She suggested that the bi- 
modality in SFHs may be intimately related to changes in 
/hi, and the transition in SF modes at ^ (2 — 3) x 
10^° Mq, found by those authors, is not a cause but an ef- 
fect of changing fni, as predicted in cold-mode accretion 
scenarios. Figure |5] illustrates these relationships by show- 
ing the averaged log f hi of galaxies which lie in different 
loci in the CMD (left panel) and the log SSFR vs. log A/* 
diagram (right panel). Contours indicate the number density 
within the q:.40-SDSS-GALEX sample in each grid point in 
the map, while the shade sc ale tra ces the mean HI fraction, 
(log /h/) . As mentioned in ^Z231 96% of ALFALFA galax- 
ies lie on the blue side of the optical-NUV CMD, whereas 
a smaller fractio n (84% ) lie on the blue side of an optical- 
only CMD (see ii4.2.3l below). Some of this difference can 
be attributed to the greater impact of shredding on the SDSS 
pipeline magnitudes relative to that of the GALEX photom- 
etry. Because of the color gradient of galaxies (outer disks 
are bluer), the shredded central redder object is identified as 
the OC. As a result, the adopted photometric object may be 
redder in (u — r) than the galaxy as a whole actually is. More- 
over, the u-band is not as sensitive and thus yields photometry 
with large uncertainties for some of the galaxies. Addition- 
ally, the {NUV — r) color is a stronger diagnostic of the SFH. 
For similar reasons, the (log fni) of grid points in regions of 
low number density should be interpreted with caution. How- 
ever the general trends (1) that the red-sequence is associated 
with low J HI and (2) that blue cloud galaxies are gas rich are 
clearly evident. At a given A/^, redder {NUV — r) colors, on 
average, indicate lower fni- Furthermore, such a variation of 
/hi along the (NUV — r) axis is more evident in the fainter 
Mr range: S{log fni) /S{NUV - r) ~ 0.75 at AU ~ -16, 
whereas S{\og fHi)/S{NUV - r) ~ 0.25 at Mr ~ -22. 
Therefore, the correlation of (NUV — r) and fni at a given 
Mr is hard to see in a sam ple with only massive galaxies, e.g. 
GASS (IWang et al.ll20l1 . 

Similarly, the right panel of Figure |5] illustrates how the HI 
fraction varies in the Af,-SSFR plane. As galaxies assemble 
their stellar mass and evolve along the star forming sequence, 
represented by the contours of high number density, their HI 
fractions decrease. In addition, for fixed A/», galaxies with 
lower SSFRs have, on average, lower fni- Again this is more 
clearly evident at the low stellar mass end, and is consistent 
with the broadening in both the /hi and SSFR distributions 
at low M^, (iHuang et al.ll2012<) . In contrast, the variation of 
(log /hi) along the log SSFR axis at a given A/» is less ev- 
ident for galaxies with > 10^ ''' M©. These two trends 
suggest that, for low A/, systems, high /hi galaxies are more 
Ukely to be starburst galaxies (defined as high SSFR galaxies). 



whereas galaxies in the high A/* regime selected by high fni 
are less likely to be starbursts. In summary. Figure |5] clearly 
demonstrates that the color, SF and gas evolution of galax- 
ies are closely related to one another, as expected. Moreover, 
the regulation of SF by Mhi is stronger in the less massive 
galaxies. 

4. THE IMPACT OF OPTICAL VERSUS HI-SELECTION 

Future surveys of HI in galaxies at intermediate to high 
redshifts that will be enabled by the next generation of 
centimeter-wavelength radio telescopes (e.g. the SKA) will 
aim to infer the gas evolution from high redshift populations 
to the local well-studied ones. It is important therefore to 
understand the fundamental properties of local Hl-selected 
galaxies, as represented by the ALFALFA catalog, and their 
biases relative to the overall galaxy population. In |3] we ex- 
amined the global properties of gas, stars and SF within the 
ALFALFA galaxies themselves. They form an Hl-rich, blue 
and less evolved population with low SFE; these characteris- 
tics are more pronounced in lower mass systems. Mhi and 
fni are linked to the SF related quantities, demonstrating the 
role that HI plays in the regulation of galaxy evolution along 
the star forming sequence. To understand how the Hl-selected 
population is biased relative to ones selected by stellar mass 
or optical flux, in this section we construct samples from both 
the a. 40 and SDSS catalogues and then compare their simi- 
larities and differences. 

4.1. Construction of control samples 

In order to ensure the galaxies contained in each of the 
optically- and Hl-selected samples are both representative of 
their respective population and fair enough to permit compar- 
ison with the other, we construct subsamples of both a. 40 
and the SDSS in the same sky volume within their common 
footprint. T he volume-limits im posed are similar to those 
discussed bv 'M artin et al.l (l2010li) . Comparable selection cri- 
teria are applied with a further requirement that acceptable 
GALEX pipeline photometry must also be available, so that 
stellar masses and star formation properties can be derived ro- 
bustly via SED-fitting. 

4.1.1. ALFALFA selected sample, Shi 

As discussed by Marti n et al.l (120 10"), radio frequency in- 
terference from the San Juan airport FAA radar transmitter at 
1350 MHz makes ALFALFA blind to HI signals in a spherical 
shell ^ 10 Mpc wide centered at a distance of ^ 215 Mpc. 
Therefore, as did those authors, we exclude 568 galaxies of 
the sample presented in fj3] which lie beyond 15000 km s^^ 
(Dsur hereafter). To maximize the overlap of contiguous sky 
coverage between the current ALFALFA and SDSS DR7, we 
consider only the two regions in the northern Galactic hemi- 
sphere (S'^OO'" < RA < 16^30"^, 4° < Dec < 16° and 
7'^40'" < RA < 16'M0'", 24° < Dec < 28°, see Figure |6};). 
Applying these restrictions yields a sample within a sky vol- 
ume of Vsur = 1.987 X 10^ Mpc^, a sky area of 1989 deg^ 
and including 7638 a.40-SDSS-GALEX galaxies. 

Next, a weight, Vsur/Vmax, is assigned to each galaxy, 
where Vmax is given by the maximum distance, Dmax, at 
which an HI source can be detected by ALFALFA, if Dmax < 
Dsur, with Vsur/Vmax = 1 for the galaxies which can be 
detected all the way outwards to the Dsur- Because the AL- 
FALFA sensitivity depends not only on the integrated HI line 
flux density, Sint [Jy km s^^], but also on the HI line pro- 
file width, 1^50 [km s"^], specifically, the fit to Sum, the 
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limiting integrated HI hne flux density that can be detected 
at S/N above 4.5 (c ode 1 and 2, 25% complete), as given in 
iHavneset all ( 120111) . is: 

, _ r 0.5 X log VKso - 1.11 - 0.202, if log Wso < 2.5; 

log asum - 1 1.0 X log W50 - 2.36 - 0.202, if log W50 ^ 2.5. 

Then, Dmax can be calculated given Sum and Mhi, based on 
the standard equation A//// = 2.356x10^ D^f^^Sint- In order 
to characterize the stellar component of the galaxies, we also 
perform a cross-match to the SDSS and GALEX databases. 
The application of such a weight scheme, or volume correc- 
tion, is equivalent to resampling the galaxies by their HI prop- 
erties (Mhi and VF50) alone, thereby reemphasizing the im- 
pact of Hl-selection. 

To further trim the sample, we drop galaxies whose weight 
Vsur /Vmax > 60, i.e., wc considcr only galaxies that could 
be detected in more than 1.67% of the survey volume. This 
cutoff corresponds approximately to a lower Mhi limit of ^ 
10*-^ M0 (see Figure|6]l. There is not a hard Mhi cutoff be- 
cause W50 also plays a role. The galaxies with Vsur /Vmax > 
60 are all relatively nearby {D„iax < 54.7 km s~^, or assum- 
ing Hubble flow, czmax < 3829 kms^^, z,nax < 0.0128), 
and are low HI mass galaxies less representative of the survey 
overall. Further motivations for applying such a weight cutoff 
include: (a) for these very local sources, distance dependent 
quantities, e.g. Mhi, have large uncertainties due to their pe- 
culiar velocities; (b) such galaxies are also underrepresented 
in the SDSS redshift sample (see below); (c) for resolved, 
patchy dwarf systems especially, the SDSS pipeline magni- 
tudes can suffer from shredding. The lowest HI mass systems 
have been considered separately in Huang et al. (2012). 

The final Hl-selected sample referred to as Shi includes 
7157 galaxies. 

4.1.2. SDSS selected sample, S opt 

To construct an optically-selected control sample out of 
the same sky volume, Vsur, we queried the SDSS DR7 in 
the same RA and Dec ranges for photometric objects which 
have valid model magnitudes and were also spectroscopic tar- 
gets. We also require them to (a) have a spectral classifica- 
tion of "galaxy"; (b) have an SDSS redshift, zsdss, deter- 
mined with high confidence; (c) lie within the same redshift 
range as Shi, czsdss < 15000 km s^^; (d) have Galactic 
extinction-corrected r-band model magnitudes brighter than 
\1 .11 . 24379 galaxies meet these criteria. Note the red- 
shifts adopted for this sample use the SDSS measurement, 
whereas that for the Shi comes from the HI line measures. 
Given the czs dss and coordinates, distances are estimated in 
the same manner as for the a. 40 sample using a local flow 
mode l for cz < 6000 km s~^, and Hubble distance other- 
wise (iHaynes e t al. 2011). Following the same procedure as 
for Shi, we searched for GALEX cross-matches, and applied 
similar SED-fitting to the UV/optical bands. 

To match the weight cutoff of the Hl-selected Shi sample, 
we also calculate weights for the SDSS-selected sample but 
here according to their optical fluxes. In this case, D„iax is the 
maximum distance at which the object, given its r-band flux, 
could be included in the SDSS main galaxy redshift sample. 
As for the Hl-selected sample, we drop galaxies with weights 
greater than 60. Given the magnitude limit of the SDSS red- 
shift survey (rrir < 17.77 mag), such a weight cut directly 
corresponds to an r-band absolute magnitude limit of ^ —16 
mag. Furthermore, since the mass-to-light ratio varies only 



mildly with color, the luminosity cut approximately translates 
to a stellar mass lower limit of ^ 10^'^ Mq. Finally, a small 
number of galaxies are removed because they are included 
in the a. 40 catalog but have been previously noted by indi- 
vidual inspection to have suspect photometry (IHaynes et aP 
I2OI 1|) . The final optically-selected sample referred to as Sopt 
includes 16817 galaxi es, o f which 34% are cross-matched to 
the a AO catalog (see m.2i . The remainder are missed by AL- 
FALFA either because they are (1) gas poor, (2) lie at a suffi- 
cient distance that their HI line flux densities falls below the 
HI sensitivity limit, or (3) for some other reason, e.g. their 
HI spectrum is contaminated by RFI or was not sampled at 
all (small gaps in ALFALFA coverage), or they correspond 
to one "child" of a shredded photometric parent object, but 
another photometric child is favored as the best match to the 
ALFALFA detection. 

We note that the distributions of weight, Vsur /Vmax, for 
both the samples highly peak at 1, and that the weight cut 
of 60 applied to each dataset is confirmed to be high enough 
to retain the main populations. Especially for the Sopt, 69% 
of the galaxies have a unit weight, i.e. can be detected out- 
wards to the edge of Vsur as we defined here. The number of 
galaxies in bins associated with a weight above 1 drops more 
rapidly among the Sopt sample, confirming that the SDSS is 
deeper than ALFALFA. 

4.2. Comparison of control samples 
4.2.1. Basic properties 

Figure |6] illustrates the comparison of quantities relevant to 
sample selection between the Shi and the Sopt populations. 
In panels (a-c), red points denote galaxies in Sopt, whereas 
blue points denote galaxies in Shi- In the histograms, red 
lines illustrate the distribution of Sopt and blue lines trace 
Shi', above each histogram, separate panels show the frac- 
tion of Sopt galaxies that are cross-matched to a. 40 in each 
bin. The numbers in each subset are indicated in panels (a) 
and (b). 

The top row contains two Spaenhauer diagrams showing, 
respectively, the r-band absolute magnitude (panel a) and HI 
mass (panel b) versus distance. The solid vertical line repre- 
sents the cz cutoff, 15000 km s"^ We use the SDSS redshifts 
to derive distances for the Sopt and ALFALFA HI v elocity for 
the Shi- As discussed also bv iMartin et al.l (1201 Oh . a survey 
must sample sufficient volume to detect very massive galaxies 
in either stellar (panel a) or HI (panel b) mass. ALFALFA for 
the first time provides a full census of Hl-bearing objects over 
a cosmologically significant volume of the local universe. 

As evident in Figure|6ja), SDSS is volume limited to A/^ ^ 
— 19 mag. The sharp lower edge of the Sopt distribution above 
£> ~' 60 Mpc results from the magnitude limit of the SDSS 
main galaxy redshift sample; as noted before, the adopted 
weight cutoff corresponds to M^ ^ —16 mag (horizontal 
dashed line). Since no limit on any optical quantity is applied 
to the Shi subset, many blue points from a. 40 show up faint- 
wards of the lower edge of Sopt, as faint as Mr ^ —14. We 
note that the blue points lying faintwards of the lower edge of 
Sopt and above D 60 Mpc are still detected by the SDSS, 
but most often only as photometric objects; hence their op- 
tical redshifts are generally unknown. To enable SED-fitting 
however, all the Shi galaxies are required to be detected in 
the SDSS; the very rare "dark" HI clouds without identified 
OCs included in a. 40 are outside the scope of this work and 
are not included in this discussion. 
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Figure |6jb) shows how Mhi increases with distance. The 
HI measures of Sopt all come from the a. 40 catalog (5653 out 
of 16817 Sopt galaxies, i.e., a cross-match rate of 34%). The 
weight cutoff applied to the Shi sample results in an approxi- 
mate limit of Mhi — 10^'^ Mq (horizontal dashed line). The 
red points below lO'*'^ M0 are still detected by a. 40, but are 
excluded from the Shi simply because o f their high weights . 
Because of the way distances are derived (iHavnes et al.l201 lb 
from the observed redshifts, an insignificant difference (0.1 
Mpc) of the median distance arises for the same galaxies in 
the Sopt-Sni overlap sample. This difference is mainly due 
to the fact that the Shi galaxies may be assigned membership 
in group whereas such information has not been applied to the 
Sopt sample. A few red outliers below the main distribution 
indicate objects without robust SDSS redshifts. 

Panel (c) shows the sky distribution of the two samples. 
Besides the large scale structure, the required availability of 
GALEX photometry also contributes to the pattern. For ex- 
ample, the Virgo region is densely covered by at least MIS 
level visits, whereas the patches of sky blank in either sam- 
ple and with regular edges arise from the lack of GALEX 
coverage in FUV and/or NUV. The distribution of distance 
is shown in panel (d). The two samples roughly coincide 
within ^100 Mpc. Above the histogram, the fraction of the 
Sopt galaxies that are cross-matched to a. 40 is shown. The 
first peak in number density coincides with the Virgo cluster 
(^•16. 5 Mpc), where Sopt clearly out numbers Shi- Since 
the two distributions agree with each other again at larger dis- 
tances ^50 Mpc, the disagreement at the Virgo distance in- 
dicates a real underdensity of gas-rich detections in Sht, re- 
flecting the well-known HI deficiency (iDavies & Lewisll 19731; 
iGiovaneUi & HavnesI 119851; ISolanes et all 120021) . Beyond 
~100 Mpc, Sopt significantly overtakes Shi, though the 
shapes of peaks or gaps still agree. This suggests that AL- 
FALFA is capable of detecting HI massive objects at large 
distances, although the survey is not as deep as SDSS. Al- 
though in a given distance bin, the least massive objects con- 
tained in Shi are even fainter than those in Sopt (i-e., the blue 
points below the lower edge of red distribution in panel a), 
ALFALFA is not as complete as SDSS at large distances. 

The distributions of Mr for both samples are shown in panel 
(e). The vertical solid line denotes the equivalent weight cut 
applied to Sopt- While Sopt peaks at a slightly fainter Mr than 
Shi, the two samples probe a similar Mr range so that their 
comparison is valid. 

The distributions of axial ratio, given by the SDSS pipeline 
measures of the exponential fit a/b in r-band, are shown 
in panel (f). Because the ALFALFA sensitivity depends 
on the HI line profile width (see Shi might be ex- 

pected to be biased ag ainst edge- o n gala xies with high a/b 
values. For example, ' West et aLl (120101) demonstrated that 
their Parkes Equatorial Survey (ES, a search through HIPASS 
cubes; Garcia- Appa doo et al.l l2009h - SDSS common sam- 
ple is slightly biased towards face-on galaxies, relative to an 
SDSS DR4 sample, with the mean log a/6 equal to 0.17 and 
0.21 for their ES-SDSS and DR4 samples respectively. How- 
ever, panel (f) shows no such obvious bias. Both Shi and 
Sopt have the same (log a/6) = 0.28. Furthermore, the cross- 
match rate even slightly rises for high a/b galaxies, with only 
a mild drop in the very last bin. Visual inspection shows that 
shredding can cause large errors in the a/b measures by the 
SDSS pipeline. The Sopt sample contains more galaxies with 
bulges making their a/b values appear to be smaller; in con- 
trast, Shi is biased against such galaxies. 



4.2.2. Internal extinction in Hl-selected galaxies 

Previously, and in many analyses of SDSS derived sam- 
ples, internal extinction is ignored. However, while the 
outer parts of galaxy disks are transparent, it is well es- 
tablished that the inner regions are optically thick at short 
wavelengths. Therefore, the neglect of internal extinction 
in disk-dominated galaxies is likely to introduce systematic 
inclination-dependent effects. In this section, we discuss (a) 
how internal extinction varies with stellar mass, (b) how in- 
ternal extinction may introduce scatter into relationships in- 
volving colors and (c) how the extinction characteristics of 
the Shi Hl-selected galaxies compare to those derived from 
an optically selected sample Sopt- In the three figures asso- 
ciated with this section. Figure |7] - 19] typical error bars on 
individual points are plotted in selected panels in as well as 
the Spearman's rank correlation coefficients. 

Estimates of internal extinction are derived from 
UV/optical SED-fitting as before. The two-component 
dust model jCharlot & Falll 120001) is incorpo rated into the 
construction of the library of model SEDs (iGallazzi et al.l 
12005ft ; the process accounts for both the diffuse interstellar 
medium (ISM) and short-lived (10 Myr) giant molecular 
clouds. Such estimates have overall larger uncertainties 
among th e red-sequence galaxie s relative to the blue cloud 
ones (e.g. ISaintonge et al.ll201 Ibh. and con t ribute to the SFR 
uncertainties. Furthermore, ISaUm et all (l2007h identified 
differences between the effective optical depth in V band, 
TV, derived from emission-line fitting and that derived from 
SED-fitting, as a function of stellar mass. Specifically, at 
lower masses, the SED-fitting-derived value is systematically 
higher than the line-fitting-derived one, but the situation 
is reversed at the high mass end. Therefore, we applied a 
Gaussian prior distribution of ry for each model, given the 
absolute magnitude and axial ratio of the individual galaxies 
(see fQ. The mean of the prior distr ibution is given by 
equation (12) in iGiovanelli et al] (11997b . which depends on 
the axial ratio and absolute magnitude, i.e. more luminous 
edge-on galaxies have larger extinctions. 

Figures|7](a-c) show plots of r-band internal extinction ver- 
sus stellar mass for the Shi galaxies. Despite the large uncer- 
tainty, internal extinction is a weakly increasing function of 
A/, in all these panels. Panel (a) shows SED-fitting-derived 
values before the prior distribution applied. A population of 
low mass red galaxies have unrealistically high Ar (SED no 
prior), because of the age-extinction degeneracy. The mean 
of the prior distribution of internal extinction, Ar (prior), is in 
panel (b). Although the Ar (prior) values of low mass galaxies 
are confined to low values, a population of massive galaxies 
have unrealistically low Ar (prior) likely due to the underes- 
timate of a/b. Visual inspection shows that shredding tends 
to describe sources as rounder in edge-on galaxies; dominant 
bulges, dust lanes and seeing effects will also lead to system- 
atic underestimates of the axial ratio. Panel (c) plots the SED- 
fitting-derived values after the prior distribution is applied, Ar 
(SED with prior). Both a lack of high mass galaxies with low 
Ar as well as of low mass ones with high Ar are evident in 
panel (c). Combining the distributions in panel (a) and (b), 
the Ar (SED with prior)-7\/, correlation is the tightest of the 
three, with a correlation coefficient rg = 0.33. Whereas ap- 
plying the prior reduces the systematic offsets of the Ar es- 
timates by SED-fitting, as well as th e SFR values, it h as lit- 
tle effect on the stellar mass values dHuang et al.ll2012l) . Fig- 
ure Ul (d) demonstrates that, as expected, the derived values 
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of Ar (SED with prior) are higher in more edge-on galaxies. 
Neglect of corrections for internal extinction will lead to the 
systematic underestimate of luminosity, so that hereafter we 
apply the SED-fitting with prior corrections. 

For the subset of the Shi galaxies (6164/7157) which are 
included in the MPA-JHU DR7 release of SDSS spectral mea- 
suremen ts Qhttp : // www, mpa-garc hing . mpg . de/S D S S/ DR74 
iBrinchmann e7al.l 120041) . we have verified that at 
M^, < the Ar inferred from the Balmer decre- 

ment and from the SED fitting using an Ar prior are in good 
agreement and the offset observed in Salimetal. (2007) is 
reduced. Above this mass the Balmer decrement leads to 
larger Ar values, but this is not unexpected as the SDSS 
spectra observe only a small region, typically towards the 
center of the galaxy where metallicities and dust attenuations 
are higher 

Another way to explore the importance of extinction cor- 
rection in a population involves examining the scatter in the 
log SSFR - color diagram, as shown in Figure[8] Results for 
individual points (unweighted) for the Sopt sample are shown 
in the left panels and for the Shi galaxies on the right, us- 
ing the colors {NUV — r)o in the top row and (u — r)o in 
the bottom, respectively. The subscript '0' in the labels in- 
dicates that the color s are corrected for internal extinction. 
As demonstrated in SSFR is an intrinsically strong 

function of (NUV — r), because NUV traces the SFR and 
the r-band luminosity is related to the stellar mass. For the 
Hl-selected population, comparison can be made directly of 
the volume-limited Shi sample shown in Figure [Ub) with 
the corresponding result, uncorrected for internal extinction, 
shown in Figuregff) for the full a.40-SDSS-GALEX. As ev- 
ident by inspection, in addition to a shift bluewards of the 
points in Figure|8jb), the dispersion about the mean relation is 
greatly reduced when the extinction correction is applied, and 
the Pearson correlation coefficient likewise improves from 
rp = —0.78 to rp = —0.86. This analysis indicates that the 
scatter in Figure |4ff) is substantially amplified by the lack of 
a correction for dust extinction. Other factors causing scatter 
include different SFHs, metallicities, as well as different pop- 
ulation synthesis models and even IMFs ( Gunawardhana et al. 
1201 ll) . However, dust extinction dominates among these fac- 
tors. Furthermore, the star-forming Shi galaxies have on av- 
erage bluer {NUV — r)o colors than the Sopt ones (median 
1.42 versus 1.74 mag). Meanwhile, the tail of red and low 
SSFR galaxies in panel (a) disappears in panel (b), again re- 
flecting the bias present in the Hl-selected population. An 
additional result of the Hl-selection is that the typical error 
bar in panel (b) is slightly smaller than in (a). 

The lower panels in Figure |8] examine the optical SDSS 
colors {u — r)o- As discussed previously in m.2\ on aver- 
age S{u ~ r)/5{NUV — r) ^ 0.6. Figure [HJc) shows that 
galaxies with SSFR < 10^^^ yr^^ have similar (u — r) col- 
ors, forming a vertical tail in these plots; this red tail is much 
less pronounced in panel (d). The adoption of (NUV — r) as 
the color index breaks down the degeneracy of (u — r) in the 
red range (u — r > 2.3) when inferences on the SFH are in- 
ferred. Furthermore, the SSFR correlates more linearly with 
(NUV — r)o than with {u — r)o, e.g. rp = —0.95 versus 
-0.87 in Sopt. 

Given a simple assumption of the dust and stellar geometry, 
a well calibrated IRX — (3 relation (IRX, infrared-excess de- 
fined as LtirI Lpuv', (3, the UV spectral slope) is sufficient 
to predict the Apuv (tightly correlated with IRX) from the 



UV color (cha racterizing the UV spectral slope, /3), in star- 
burst galaxies jCalzetti et al .1119941) . However, the IRX - (3 
relation in normal star-forming galaxies has a shallower slope 
with larger scatter (Salim et al. 200 71) - This r esult may be due 
to differences in the SFH (Kon g et an i2004'). or dust geome- 
try ( Cor tese et al. 2006). Therefore, the loci of galaxies on an 
Apuv versus {FUV — NUV) plot give an indication of the 
dust extinction behavior. 

Following this approach, the two samples Shi and Sopt are 
compared in Figure|9] with results for Sopt on the left and Shi 
on the rig ht, b oth before and after applying the weight correc- 
tion (see ^4Ai in the upper and bottom rows respectively. The 
red dashe d line corresponds t o the fit to star-forming galaxies 
derived in lSaUm et al.l ( 120071) . based on a typical local SDSS- 
GALEX cross-matched catalog. Despite the large scatter, it 
is on average in close agreement with the distribution of the 
Sopt galaxies. Galaxies closer to the fit have overall higher 
weights so that the correlation appears to be slightly tighter 
in panel (c) than in (a). However, the distribution of the Shi 
galaxies is offset from the fit in panels (b) and (d), toward the 
lower Apuv side, i.e. for a given {FUV — NUV) color, the 
Hl-selected galaxies on average have lower extinctions. This 
result is also confirmed by the generally lower Apuv of the 
Shi sample overall, with a median value of 1.31 mag, rela- 
tive to that of Sopt, 1.4 6 mag. Such a dev iation of the main 
trend from the SF-fit of lSalim et"an (l2007h suggests that the 
Hl-selected galaxies have different SFHs or dust geometries. 

Unfortunately, the metallicity of the stellar population is 
poorly constrained by the SED-fitting. Nevert heless, given 
the co rrelation between dust and metallicity (e.g. lDraine et al.l 
l2007h . the lower extinction infers that the Shi galaxies 
have lower metallicity. Besides t he well-known mass- 
metallicity (gas - phase) relation (e.g. iTremonti et al.l |2004), 
iMannucci et al.l (120 10 *) demonstrate that at lower A/»(< 
IO^^'^Mq), metallicity decreases sharply with increasing 
SFR, while at high stellar mass, metallicity does not depend 
on the SFR. Given the b ias tow ards high SFR of the Shi sam- 
ple at a given (see ! j4.2.4l i. the bias towards low metallic- 
ity and low extinction is expected. The gas-phase metallic- 
ity measures (oxygen abundance) from the MPA-JHU DR7 
release (Tremonti et al. 2004) are available for 4211/7157 of 
the Shi galaxies and 10311/16817 of the Sopt galaxies. In 
addition to the caveat of the small SDSS fiber aperture, the 
requirement of being an SDSS spectroscopic target may re- 
duce the difference between the two subsets. The mean 
12 + \og{0/II) is only slightly higher among the Sopt galax- 
ies (8.74) than in the Shi galaxies (8.71). The overall lower 
extinction among the Shi galaxies is consistent with the 
idea that Hl-selected galaxies are relatively gas rich and less 
evolved, with more likely gas infall, lower SFE and metallic- 
ity, less dust and thus lower extinction. 

4.2.3. Distribution in the intrinsic CMD 

The bimodal distribution in the optical CMD dBaldrv et al.l 
120041) has been interpreted as an evolutionary sequence, with 
the blue-cloud galaxies growing through mergers and the con- 
sumption of gas and later migrating to the red sequence. This 
evolutionary scenario also predicts that galaxies evolve from 
a state of low stellar mass, hig h SS FR and high HI fraction to 
the opposite. As discussed in ^3A\ the HI fraction is lower in 
redder galaxies leading to the result that an Hl-selected sam- 
ple like Shi is biased against red galaxies whereas those are 
commonly included in, and may even dominate, an optically- 
selected one like Sopt- Here, we use the CMD to quantify this 
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bias. 

Fi gure fTOl shows the CMDs of the Sopt (left panels) and the 
Shi (right panels); before and after applying weight correc- 
tions in the upper and bottom rows respectively. In these plots, 
absolute magnitudes (Af, )o and (u — r)o colors are corrected 
for internal extinction as given by the SED-fitting with the 
prior applied. A similar CMD for the whole a.40-SDSS pop- 
ulation but without corrections can be found in Figure 7 of 
iHavnes et al.l (1201 ll) . The dash-dotted curve is based on the 
best fit to the division of the red and blue populations derived 
by Baldry et al. ( 2004) but with shifts toward bluer colors and 
brighter Mr appUed in the amount of the mean extinction of 
the Sopt sample. 

As is obvious in panels (a) and (b), the red sequence is far 
more pronounced in Sopt than in Shi- Among the Sopt galax- 
ies, only 68% lie on the blue side of the division, whereas the 
percentage of blue galaxies in the Hl-selected Shi popula- 
tion is as high as 84%. The small population of 'red' galax- 
ies with HI represents candidate objects which (i) have re- 
cently migrated onto the red sequence retaining some resid- 
ual gas, or (ii) will transit back to the blue cloud via late gas 
re-accretion (Cortese & Hughesl l2009h . Note that most of the 
red Hl-bearing galaxies optically luminous. The HI mass of 
the faint red galaxies is usually so low that their HI line flux 
densities are below the sensitivity limit of ALFALFA. 

In the Vsur/Vmax Weight corrected diagrams (panels c and 
d), the peaks of number density shift to fainter {Mr)o in 
both samples. Galaxies densely occupy the faintest (Afr)o 
bins, especially in the Sopt sample. This result is consistent 
with the optical lum inosity function of the blue population 
dBaldrv et al.ll2004h . A peak of the number density on the 
red sequence around (Af,.)o ^ —21 also coinc ides with the 
maxim a of the luminosity function derived by iBaldrv et al.l 
(120041) . However, a second peak appears on the red sequence 
at the faint end. Note that the data points in Baldry et aT| 
(|2004) also suggest a rise in number density in the faintest 
bins of the red sequence luminosity function (see their Fig- 
ure 7), though those authors ignore these objects when they 
fit the Schechter function. In agreement with this. lHogg et al.l 
(|2b03) reported a non-mono tonic density trend along the red 
sequence, seen as a dip in the typical density for intermediate- 
mass red-sequence galaxies despite higher densities at higher 
and lower masses. Similar to the finding of a sudden broad- 
ening in the SS FR distribution faintwards of Mb ^ —15 
dLee et al.ll2007h . a sudden broadening in the color distribu- 
tion faintwards of {Mr)o ^ —16 is seen in the Sopt sample. 
This increased spread suggests that, unlike the galaxies of in- 
termediate mass in which the SF is mainly regulated by stellar 
mass and for which the blue population dominates in number, 
the gas in dwarfs can be easily removed so that SF suddenly 
quenches driving their migration onto the red sequence. 

Figure [TOf c) shows a CMD for an optically-selected pop- 
ulation similar to Figure 2 of iBaldry et al.i (l2004l) . How- 
ever, Sopt is limited to cz < 15(3 km whereas 
the low-redshift sample iBaldrv et al.l (120041) studied has 
1200 kms-i < cz < 24000 km s'^ The imposition of 
the volume cut here clearly decreases the presence of lumi- 
nous red galaxies; the local volume is dominated by blue star- 
forming galaxies. This effect is also enhanced by the addi- 
tional UV selection applied to the Sopt sample. 

The near-absence of the red sequence among the Shi popu- 
lation is clear in panel (d); as we have noted before, a blind HI 
survey like ALFALFA is highly biased against red sequence 
galaxies. Additionally, in both populations galaxies appear 



to be bluer as their Mr gets fainter, though with large scatter 
Such a trend is weaker but still visible after internal extinction 
corrections are applied. Therefore, the slope is not only due 
to extinction, but is intrinsic. 

Besides comparing the CMDs of Sopt and Shi, we can also 
study the impact of HI selection through an examination of 
the fraction of Sopt galaxies that are cross matched to a. 40. 
As was discussed in ^^A] 34% of the Sopt galaxies are cross- 
matched to the a. 40 catalog. However, we note that such a 
fraction is a lower limit of the Hl-detection rate of optical-UV 
selected galaxies in a. 40, because some Sopt objects may be 
a shredded photometric object of a gas bearer while another 
piece is cross-matched to the a. 40 entry. The fraction of thes e 
objects should however not be large (see IHavnes et al.|[201 Ih . 

Figure [TT] explores the fraction of the Sopt population in- 
cluded in the a.40 within the CMD (left) and the SSFR- 
M^ diagram (right). In the left panel, the NUV-to-optical 
CMD with internal extinction corrections is color coded by 
the fraction of galaxies in the optically-selected Sopt sample 
which are also included in the q.40; the contours trace the 
total number density of the Sopt galaxies. Clearly, an HI sur- 
vey like ALFALFA is more efficient at detecting blue galax- 
ies than red ones. Especially at very bright end of the blue 
cloud (lower left corner), the detection rate is close to 100%. 
Furthermore, the ALFALFA population also detects the ma- 
jority of optically faint blue galaxies with the detection rate 
again approaching ^100% (lower right corner). The latter 
have the highest HI fractions (see also Figure [S]). HI survey 
is most efficient in detecting galaxies with large diameter and 
high Mhi /L values ( Garcia- Appadoo e t al. 2009 ). The gen- 
eral shape of the color scale variation above a cross-match 
rate of ^40% resembles that of the number density variation 
on the bluer side of the blue cloud. In contrast, starting from 
the red edge of the blue cloud and throughout the whole red 
sequence region, the rate of detection by ALFALFA drops be- 
low 10% and even to ^0%. These results again demonstrate 
that Hl-selection is highly biased against galaxies on the red 
sequence and furthermore, that such bias begins to become 
pronounced on the redder portions of what still is considered 
the blue cloud. The precise location of this limit to some ex- 
tent reflects the 'fast' nature of ALFALFA (only 40 seconds 
integration per beam); a deeper survey of comparable solid 
angle would probe to lower gas fractions and hence deliver a 
higher cross-match rate. 

4.2.4. Star-formation sequence 

The scaling relation that the SFR increases with among 
the a.40 galaxies is demonstrated in i j3.2l (Figurel4h). Figure 
Ud) similarly shows the well confined star forming sequence 
in flie SSFR-M ^ diagram. Based on the ES-SDSS sample, 
IWest et al] (120091) suggested that the significant offset towards 
bluer colors, induced by Hl-selection can be explained by en- 
hanced recent bursts of star formation. Furthermore, rising 
SFHs, i.e. high 6— parameter, are required to explain the col- 
ors of Hl-selected galaxies bluer than g — r < 0.3, which may 
result from gas infall and enhanced SF subsequently. Here, 
we compare the distributions of the 5*^/ and Sopt galaxies, 
on and off the star-forming sequences to explore further the 
impact of Hl-selection. 

Figure [12] examines the SFR versus A/* scaling relation 
for the optical (left) and HI (right) selected samples. The up- 
per row shows the individual galaxies while the lower traces 
the Vsur/Vmax Weighted distributions. The cyan dashed line 
at M* > 10^° Mq shows the fit to this relation obtained by 
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ISchiminovich et all (120101) . based on the high Af* GASS sam- 
ple in all four panels: 

log(5Fi?(M,)) = 0.15 log - 1.5, logAf* > 10.0. 

Like Sopt, GASS is essentially an SDSS selected sample but 
only includes the high stellar mass end. The blue diamonds 
and lines represent \og{SFR) derived for Sopt or Shi in each 
stellar mass bin. Note these are different from the average 
trends in Figure|2]or Figure|4] The values are weighted aver- 
ages in the weighted panels. The numbers of galaxies in each 
stellar mass bins are listed at the top of panels (a) and (b). 
Typical error bars are shown in panels (a, b). 

The main distribution, consisting of the star forming galax- 
ies, is associated with the parameter space inside the con- 
tours. In the un-weighted diagrams, the GASS fit is consistent 
with Sopt in the high mass range (panel a), but falls below 
the average of Shi (panel b). This offset indicates that even 
when only the star forming sequence is considered, an Hl- 
selected population has a higher SFR overall, because Shi is 
not complete even within the blue population (see also Figure 
fTTT i. Similarly, the blue diamonds in the weighted Sopt plot 
roughly coincide with the GASS fit (panel c). At the same 
time, the \og{SFR) of Shi galaxies becomes systematically 
larger than the GASS fit with increasing Af, (panel d). 

Below the locus of star forming galaxies, the Sopt sample 
contains a large population of passive galaxies visible as the 
cloud of points extending to very low SFRs (panel a). How- 
ever, this population is barely seen in the same diagram for the 
Shi galaxies (panel b), i.e. an Hl-selected sample is highly 
biased against non-star-forming galaxies. 

Because Sopt is truncated at the l ow m ass end due to the ap- 
plied Vsur/Vmax Weight limit (see ^^A^, Shi probes to lower 
A/, and thus lower SFRs on average. Despite this dwarf tail, 
the average SFR is still slightly higher in Shi- The median 
SFR is 0.4 M© yr"i in Sopt versus 0.6 Mq yr^i in Shi- 
In the weighted plots, the correlation extending to the low A/» 
range is better illustrated (panel d). There is a hint of a change 
of slope in the SFR-M^, scaling relation below the SF tran- 
sition mass Af, ~ lO^'^Af© noted earlier (! j3.2l and equation 
[8]l, with the slope at high Af, being shallower, especially evi- 
dent in the Sopt plots. We note that the trend appears to flatten 
again in the lowest A/, bin in panel (c) because of the weight 
cut. To be more specific, at a given Mr, the bluer galaxies 
with higher SSFRs have on average lower stellar masses, ac- 
cording to the mass-to-light ratio versus color relation. There- 
fore, the uniform A/^ cutoff applied on the Sopt sample makes 
the distribution of Af, extend to lower limits in higher SSFR 
bins; this effect is visible in the Sopt diagrams (Figure [T2h . 
c; see also the right panel of Figure fTTT i. This effect artifi- 
cially results in the higher \og{SFR) value in the lowest stel- 
lar mass bin. Furthermore, the lowest Af, bin only contains 
70 galaxies, so that the average is even less reliable for sta- 
tistical reasons. Such a sudden flattening is not seen within 
the Shi population (panel d), which probes to lower A/, and 
with better sampling. 

Figure [13] shows similar plots of the SSFR - Af, correla- 
tion. The cyan dashed line again denotes the \og{SSFR) val- 
ues tabulated in Schiminovich et al. (2010), derived from the 
high A/, GASS sample. The contoured re gion roughly repre - 
sents the star forming sequence defined in lSaUm et al.l(l2007h . 
The red dashed line is the fit of the sequence to the blue galax- 
ies (NUV — r < 4) only as derived by those authors for a 
typical local SDSS-GALEX cross-matched catalog, with the 



majority falling in the stellar mass range of 10 — 10 M©: 

1 QQru _ /-0.17(logA/, - 10) - 9.65, logAf, < 9.4 
iog:3:3fn-~ \ -0.53(logA/, - 10) - 9.87, log Af, > 9.4. 

For comparison, the blue diamonds in Figure [13] are derived 
for our datasets but also considering only the blue galaxies 
(NUV — r < 4). Note these fits are different from the one 
given as equation [8] shown in Figure |4jd). 

Similar to what was seen in Figure [12] the Shi popula- 
tion has on average higher SSFRs than the Sopt galaxies. The 
median SSFR is lO-^-^^ yp-i for Sopt versus lO"^-^^ yj.-i 
for Shi- The breakdown of the star forming sequence above 
stellar mass ^ 2 x 10^° M© is only evident among the Sopt 
population, whereas comparably low SSFR galaxies are al- 
most entirely absent from the Shi sample. The GASS result 
agrees with the contours of Sopt above 10^° M© (panel a and 
c), but lies below the average of Shi (panel b and d). 

In panel (c), the contours, associated with the high num- 
ber density regions, that fall far below the main star forming 
sequence (log SSFR < —11) mostly show up only at the 
highest and lowest Af* ends, suggesting that the scatter in the 
SSFR distribution is larger at both A/* ends than in interme- 
diate mass range. However, the effective Mr cutoff applied 
to Sopt inhibits us from a more convincing conclusion on the 
possible breakdown of the star forming sequence at the low 
mass end. 

At the high A/* end, the averaged SSFRs are systematically 
higher in the Shi population than in the Sopt one. However, 
when only blue galaxies are considered in either sample, the 
discrepancy between the two is reduced so that the trends in- 
dicated by the blue diamonds are comparable in both. The 
blue diamonds are in good a greement w i th the definition of 
the star forming sequence in lSalim et alJ (1200 7). particularly 
for the high mass galaxies (Af* > 10^-^ M©) in the weighted 
plots, whereas the Shi distribution is slightly offset to higher 
SSFRs in an un-weighted plot (panel b). In ad dition, note 
that w e obtain a steeper slope than the blue fit in ISalim et al] 
( 12007b for low mass galaxies, even in the Sopt sample. Again 
because of the cutoff applied to Sopt, the averaged SSFR value 
in the lowest mass bin is artificially high. 

Finally, we inspect the cross-match rate of Sopt to a. 40 over 
the SSFR - Af* diagram in the right panel of Figure [TT] The 
detection rate is the highest among galaxies with high SSFRs 
at both high and low AI^, ends (above ^ 60%) and is close 
to the overall average (^ 40%) throughout the high number 
density region along the star forming sequence. However, it 
drops to below ^20% from the lower edge of the star-forming 
sequence, and the slope of the low match rate division, which 
corresponds to a cross-match rate of ~ 25% (yellow band 
in this plot), coincides well with that of the star-forming se- 
quence. The SF properties of a galaxy have a significant im- 
pact on its likelihood of detection by an HI survey. 

We again note that the blue diamonds in Figure \l3\ h) are 
slightly above the red dot-dashed line suggesting that not 
all star-forming galaxies are detected by ALFALFA. From 
the fact that the HI detection rate of the galaxies in the SF 
sequence is far from 100% (see right panel of Figure fTTI) . 
one may naively conclude that not all star-forming galaxies 
contain HI. However, this conclusion is not necessarily cor- 
rect. For example, given the beam size of the ALFA frontend 
(^3.5'), confusion can contribute to a lower effective detec- 
tion rate, e.g. when a close pair of star-forming HI bearing 
galaxies is unresolved. In fact, the red dot-dashed line is in 
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excellent agreement with the blue diamonds in Figure [T3l d). 
indicating that the mean trend of the star forming sequence 
can be well reproduced by the a. 40 galaxies after volume cor- 
rection. An Hl-selected sample can effectively trace the aver- 
age of the star forming sequence despite the non-detection of 
some star-forming galaxies by ALFALFA (which is by de- 
sign a short-integration time survey). In order to predict the 
SSFR of an Hl-selected sample, one should still use equa- 
tion (O rather than the average scaling relations derived in 
this section. The latter relations, as well as the star forming 
sequence derived in Salim et al. (2007) rely on additional in- 
formation, e.g., UV magnitudes or emission lines, to define 
the star-forming galaxies. However, the {NUV — r) < 4 (or 
similarly log SSFR < —11 yr^^) is crude, since the SSFR 
generally decreases with increasing Af*. Furthermore, there 
are actually some very low SSFR galaxies detected in HI (e.g. 
iHalien beck et al. 2012). Direct adoption of the SF sequence 
defined in , Salim et al.. C2007.) will overpredict the SSFR of an 
Hl-selected sample. 

5. THE HOST HALOS OF HI-SELECTED GALAXIES 

The results of the previous sections suggest that HI blind 
surveys detect in abundance the star-forming population but 
are highly biased against red sequence galaxies. Conversely, 
nearly all star-forming galaxies contain HI. Because AL- 
FALFA is both wide area and sensitive, it samples the HI 
mass function over 4 orders of magnitude 7 < log Mhi < 11 
with a mean of 9.56 dex (Martin et al. 2010) and, while not as 
deep as the SDSS, both surveys probe the same range of Mr- 
The brightest and reddest galaxies are missing from a. 40, but 
an Hl-selected sample provides an important perspective on 
the star forming sequence. While galaxies in rich clusters are 
well known to be HI deficient because of the secular evolu- 
tionary processes (ram pressure stripping, thermal viscosity, 
harassment), Hl-bearing galaxies dominate the extragalactic 
population in the low density field. As discussed in ^3.2\ the 
ALFALFA galaxies have, on average, lower SFEs, and equiv- 
alently therefore, longer tji than the optically selected popu- 
lation. Both this work and the study of the lowest HI mass 
dwarf population by (Huang et al. 2012) suggest that overall, 
the Hl-selected population is less evolved and has higher SFR 
and SSFR but lower SFE and extinction than one selected by 
optical brightness or stellar mass. 

Perhaps the most surprising result of the ALFALFA sur- 
vey has been the detection of a large number of galaxies with 
very high HI masses, log Mhi > 10.3, including a signifi- 
cant number with abnormally high HI gas fractions (for their 
stellar mass), dubbed the High HI Mass (HlghMass) galax- 
ies. The existence of this population, albeit rare, begs the 
question: how can such massive disks retain their huge gas 
reservoirs without forming stars? One explanation is that, 
while recent inflow of HI gas has taken place, SF in that 
gas has not yet been triggered. Alternatively, the overall low 
SFEs characteristic of the ALFALFA population may be ex- 
plained by the semi - analy tic models of galaxy formation in 
iBoissier & PrantzosI (12000 *). which predict that more extended 
disks with larger scale length and lower stellar surface den- 
sity are associated with dark matter halos of high spin pa- 
rameter A. Selected examples of very massive but gas-rich 
galaxies have been studied in recent years (e.g.. lPortas et al.l 
l2010i) . with most belonging to the extreme category of low 
surface brig htness (LSB) galaxies k nown as the "crouch- 
ing giants" (Disnev & Phillipps' 1987) or "Malin 1 cousins" 
dBothun et al...l987;.lmpev & Bothun.1997.) . Their character- 



istic huge size, low optical surface brightness, star formation 
rate and low metal abundances may be explained by the large 
angular momentum and thus the low surface density of their 
ISM; star formation in such disks is suppressed according to 
canonical star formation laws. However, we note that many 
of the HlghMass galaxies are not crouching giants. Instead of 
being quiescent objects like Malin 1, most of the HlghMass 
galaxies have blue disks and exhibit healthy on-going SF." 

The discussion throughout this paper is based on the sce- 
nario of hierarchical galaxy formation through mergers. Al- 
ternatively, cold mode accretion (Keres et al. 2005) can build 
up gaseous galaxies rapidly at high redshift (z > 1). In the 
local universe volume probed by ALFALFA, it is likely that a 
large fraction of the gas accretion happens in a slower "hid- 
den" mode, e.g., from the ionized hot corona or driven by the 
galactic-fountain process (Marinacci et al. 2010). In a future 
work, we will investigate the angular momentum distribution 
in the HlghMass galaxies to investigate whether the major- 
ity of the HI in local gas-rich galaxies reaches the assembling 
halo through filaments of cold flow at high redshift or results 
from gradual cooling out of the hot coron a. 

According to the Boissier & PrantzosI (1200 0') model, for a 
given halo mass, galaxies whose halos are characterized by 
different values of A have similar stellar masses and current 
SFRs. However, galaxies in high A halos exhibit higher gas 
fractions, lower metallicities and bluer colors than those in 
low A halos. In this picture, the halo mass controls current 
absolute values while the spin parameter determines mainly 
the shape of the SFH. Blue colors indicate high SSFRs, i.e., 
suppressed SF in the past relative to the current SFR. Massive 
compact disks have the shortest time scales of gas infall, rapid 
early SF, and thus are dominated by old stars today. In con- 
trast, the low mass galaxies reside in high A halos can even 
have a rising SFH instead of an exponentially decaying one; 
their timescales of gas infall, as well as SF, are long. This pic- 
t ure a grees with the downsizing scenario of galaxy formation 

(E3. 

Meanwhile, the low stellar surface density found in gas-rich 
systems results in a weaker gravitational field, and, by exten- 
sion, a lower mid plane g as pressure and a higher fraction of 
dif fuse HI g as (Ostriker et al. 2010). The model constructed 
by iFu et alT (2010). tracking both the atomic and molecular 
gas in disk galaxies, predicts a low HI to H2 conversion effi- 
ciency in the high A galaxies. The HI, stellar and SF proper- 
ties of the q:.40-SDSS-GALEX galaxies, relative to the other 
local samples, therefore suggest that the Hl-selected galaxies 
are biased towards ones in high A halos, and thus are blue, in- 
efficient in SF and less evolved. Nevertheless, they have com- 
parable or even slightly higher current SFRs relat ive to the 
ones in the low A halos with similar halo mass (see !j4.2.4l i. In 
this section, we test this hypothesis. 

5.1. Spin parameters derived from the Tully-Fisher relation 

In their analysis, Boissier & Prantzos (2000) adopted the 
scaling properties derived by Mo et al, (1998J, in the frame- 
work of the CDM scenario for galaxy formation. In this sce- 
nario, primordial density fluctuations give rise to dark matter 
halos of maximum rotational velocity Vhaio, with the mass of 
the halo Mhaio oc Vj^^i^. Inspired by theoretical studies to 
break the degeneracy, a second paramete r, the spin param eter 
A , is introduced to describe the halo (e.g. lMo et alJl998h . de- 
fined as A = ,/|i?|^/'^G'^^Af^^^^'^, where J is the angularmo- 
mentum and E is the total energy of the halo. Observationally, 
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A is seen to determine p roperties such as color, dis k thickness 
and bulge-to-disk ratio dHernandez & Cervantes-S odi 2006). 
Within the halo, baryonic gas condenses later and forms the 
stellar disk with mass and characterized by an exponen- 
tial surface density profile with disk scale length R^. Those 
authors also assumed that the density profile of the dark mat- 
ter halo is isothermal and responsible for establishing the flat 
disk rotation curve Vrot- Under the further assumptions that 
the potential energy of the galaxy is dominated by that of 
the halo and that it is a virialized gravitational structure, etc., 
iBoissier & Prantzos (2000) related the quantities describing 
the halo to those describing the disk, and express the spin pa- 
rameter as 



A 



GMhalo 



See also [Herna ndez & Cervantes-SodH (120061) for similar 
derivation of equation (|9|l. The only unobservable quantity 
Mhaio can be replaced by Mbaryon/ F where F is the bary- 
onic fraction, or, alternatively, M^,/Fd where Fd is the stellar 
disk mass fraction. 

Based on this framework [Hernandez et al.l ('2007) analyzed 
the empirical A distribution of samples taken from the SDSS. 
Because their sample lacked direct measures of Vrot, those 
authors invoked the Tully-Fisher relation to infer Vrot from 
the disk luminosity. Furthermore, a disk mass Tully-Fisher 
relation and a constant disk mass fraction Fd — 0.04 are 
assumed to eliminate the Mhaio te rm. As a first approach 
for direct comparison with Hernan dez et al.l ((2007), we repeat 
the method used in that work to derive the distributions, sepa- 
rately, of the A spin parameter for the the Hl-selected Shi and 
the optically-selected Sopt galaxies. We adopt the A estimator 
proposed bv .Her nandez et al. (2007]), 



A = 21.8 



(V^™t[kms-i])3/2 



(10) 



[Hernandez et alj (I2007h adopted the i?-band Tully-Fisher re- 
lation derived bv lBarton et al.l (12001 ') to assign V-ot to the ob- 
served galaxies; since this relation is valid only over the ab- 
solute magnitude range —20 > Alfj > —22.5, they restrict 
their analysis to that range. To minimize errors due to internal 
extinction, they trim the sample to leave only spiral galaxies 
( lPark& Choill2005l) having b/a > 0.6. 

The i?-band luminosity is inferred from the SDSS bands 
based on Lupton (2005), 



(M, 



R)Q 



ro - 0.2936 x (r - i)o - 0.1439. 



We applied the same absolute magnitude cut —20 > 
(Mfl)o > —22.5 to both the Shi and the Sopt samples. 
Within our distance range, the SDSS is volume limited to 
Mr ^ — 19 mag (Figure|6^). We also require that the r-band 
light profile should be better fitted by an exponential profile 
than a deVaucouleurs profile to b e sure we were i ncluding 
mainly disk galaxies. Since Herna ndez et al.| ( |2007|) ignored 
internal extinction but we have corrected for it using the SED 
fitting estimate of An, the b/a requirement in our case is less 
strict, b/a > 0.35. Such a cut is adopted so t hat the 6/o distri- 
bution is close to being flat above 0.35. Like Hernandez et"an 
(12007.) . we convert the SDSS r-band axial ratio to inclination 



cos^ i 



{b/af 



adopting an intrinsic axial ra tio of the disk of go = 0.18 as 
proposed bv ICourteaul ( Il997h . We note, however, that since 
that study was restricted to a relatively small number of Sb- 
Sc galaxies, a more conventional value of go = 0.2 may be 
more appropriate (see i j5.2l ). 

Firs t, we confirm that the TF relation [Hernandez et alJ 
(120071) used is systematically consistent with the Shi galax- 
ies, albeit with significant scatter. The HI line width 
of the Shi galaxies, W50, is converted to Vrot through 
Vrot — (W^5o/2/ sini), where the small broadening effect 
(^ 5 km s^^) of turbulence and non-circular motions on the 
observed HI linewidths is ignored. The average trend evident 
between (Af^Jo and Vrot agrees with the relation, i.e., ap- 
plying such a TF relation should have little effect on the mean 
value of the A distribution. The difference of the mean A value 
derived here as compared with that derived next in i i5.2l will be 
mainly due to the different assumptions of the baryon fraction 
F. 

Next, we obtain the A distributions of the Shi and the 
S„r,t samples, accordi ng to equation ( fTol i. and following 
[Hernandez et alJ (120071) . we confirm that they both are well 
fit by a lognormal function: 



P(Ao,CTA;A)dA = 



1 



exp 



ln'(A/Ao) 



2<yl 



dA 



The best fit parameters are Aq = 0.0525, ax — 0.422 for 
the Shi distribution of the HI selected galaxies Shi and 
Ao = 0.0489, CTA = 0.446 for the SDSS-selected Sopt distri- 
bution. Because we correct Mr for extinction, our An valu e 
is slightly below the value found bv [Hernandez et al] (l2007h . 
Ao = 0.0585. While suggestive, this difference of Ao between 
the two samples Shi and Sopt, under the assumption of a con- 
stant disk-to-total mass fraction Fd, is small. The various cuts 
and selection effects applied to both samples, e.g. detected in 
UV, better fit by an exponential profile, small dynamic range 
of Mr, may have already minimized the possible difference 
in Ao within the full sample. Given that the remaining galax- 
ies in the Shi sample are brighter on average than those in 
the Sopt sample, i.e., Vrot is higher on average, the slightly 
higher Ao found by this analysis is due to the fact that Shi 
galaxies have on average larger Rd values (see equation [TOb. 
HI selection instills a bias towards more extended galaxies. 

5.2. A distribution of the parent population 

The SDSS galaxies analyzed by [Hernandez et alJ ([2007[) 
lack direct measurements of Vrot, whereas the q;.40-SDSS- 
GALEX sample has homogeneous W50 measurements for all 
its galaxies. We hence improve the derivation of the A distri- 
bution for all the a.40-SDSS-GALEX galaxies in this sec- 
tion. Assuming that an Hl-selected sample has little contam- 
ination from galaxies which are not disk dominated, we drop 
the requirement of a higher probability of an exponential fit 
and adopt an intrinsic go = 0.2. Because HI line widths are 
available, we no longer rely on the Tully-Fisher relation to in- 
fer Vrot, and for that reason, the {Mr)o limit is also dropped. 
Most importantly, we estimate A directly from equation (|9|l 
assuming a non-constant F to derive Mhaio- 

Following Hernandez et al. (2007), the preceding assumes 
that all halos are associated with the same disk mass frac- 
tion Fd = 0.04. However, based on abundance matching 
between the observed stellar mass function and dark mat- 
ter halo mass function derived from CDM simulation s, it is 
quite well established that Fd is not a constant (Behroozi et alj 
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120 lOl and references therein). Instead, the distribution of 
as a function of Mh aio peaks around Mhaio ^ 10^^ 
(e.g. ' Guo et alj l20Tol) . i.e., low mass galaxies, correspond- 
ing to M* < 10^°'^ Mq, retain fewer baryons as their halo 
mass declines. The trend turns over abo ve M^. ~ l O^*^-^ Mq. 
Based on the a. 40 catalog, Papastergis et al.l (l20TT ) have stud- 
ied how galaxies with different Vrot populate dark matter ha- 
los, by applying similar abundance matching to the velocity 
width function, under the assumption that rotational velocity 
well traces the dark matter mass. We adopt the Vhaio-Vrot 
relation from Table 1 of Papastergis et al. (2011), assuming 
Vhaio,r,iax = 360 kms"^ The resulting Vhaio is converted 
int o the virial mass of the halo following the simulation result 
of lKlvpinetal.l (l20Tlh: Vuaio = 2.8 x l{)-\Mhaiohf-^^'' . 
For the disk scale length Rd, we use that measured by the 
SDSS pipeline in the r-band. 

Among the 9417 a.40-SDSS-GALEX galaxies, 1829 are 
dropped because of the adopted axial ratio cut {b/a > 0.35 as 
in i jS.lb . and a further 130 are dropped because they lie outside 
of the valid range of the Vrot-Vhaio relation (16 km s^^ < 
Vrot < 431 km s~^). In the end, A values are estimated for 
7458 a. 40 galaxies; their distribution is shown in the Fig- 
ure [T4fa). The black solid line traces the normalized PDF 
for the a. 40 galaxies. The red dash-dotted line in Figure 
[T4T a) shows the best log-normal fit with Aq — 0.0585 and 
a\ = 0.446 from lHernandez et all ( 120071) and the blue dashed 
line is the bes t fit to the a. 40 galaxies with Aq = 0.0929 and 
ax = 0.875. iHernandez et al.l (I2007h noted that values they 
found for both Ao and ax are in good agreement with results 
from cosmological simulations. However, the A distribution 
of the a.40-SDSS-GALEX galaxies has a mean well above 
this value, A — 0.0852. Furthermore, the distribution for 
the a.40-SDSS-GALEX sample has a much wider disper- 
sion than the previous result, arising mainly from the fact that 
we adopt different F values according to the Vrot, rather than 
treating it as a constant. As a result, the distribution of spin 
parameters of the Hl-selected population is no longer well fit 
by a lognormal function. The KS statistic implies that the 
probability that the A distribution for the Hl-selected popu- 
lation is drawn from the same underlying distribution with a 
PDF of the best-fit lognormal function with Aq ^ 0.0929 is 
negligible. 

The technique of assigning Mhaio is valid only in a statis- 
tical sense, rather than as a concrete measure for individual 
galaxies. For this reason, we look only for a general trend in 
the dependence of the luminous component on the A of the 
halos. Figure fl^l b) shows a fni-^I* diagram similar to Fig- 
ure |2c), but color coded by the mean A value of galaxies in 
each grid. The black contours outline the distribution of the 
7458 a. 40 galaxies in the high number density region. A trend 
is evident that, in a given A/, bin, the gas rich galaxies with 
higher fni on average reside in halos with higher A. Mean- 
while, along lines of constant Mhi (straight lines with slope 
of ' — 1' in the Jhi-^I* diagram), (A) barely changes with 
increasing Af,. Taking into account the mean scaling relation 
defined by the a. 40 galaxies (blue diamonds in Figure|2};) rel- 
ative to the similar relations obtained for other samples, this 
result reinforces the hypothesis that the Hl-selected galaxies 
favor high A halos. The fni-^I* relations confined by the 
'Hl-normal' or the 'outside Virgo' galaxies in Figure|2jc) are 
close to the well accepted A ^ 0.05 region, the yellow band 
in Figure [T4T b). 

We note that such a clear trend is only weakly due to the 
overall larger Rd, i.e., lower surface brightness, which the 



more gas-rich galaxies have in a given A/» bin, but is in- 
stead largely due to the rapidly increasing F — {Mhi + 
M^)/ Mhaio with Mhi in a given bin. However, it is 
known that the increase in F with Mhaio is reversed near 
M, 10^°-^ M0. In fact, the variation of (A) along the 
/hi axis disappears around A/, ^ lO^*''^ Mq in Figure 
fT4l" b'). In addition, we note that within the luminosity range 
of -20 > {Mji)o > -22.5, the stellar disk fraction, Fd, 
inferred by this method spreads over a range with a median 
of 0.04; the Mhaio-F relation crosses the turnover around 
{AIfi)o ^ —21.5. Despite the systematic variation of Fd we 
ad opt, t he median is consistent with the assumption followed 
infe] 

Because the Vrot-Vhaio relation we have used is less con- 
strained at the highest masses, it barely reproduces the trend 
of decreasing F with Mhaio above A/* IQ^^-^ Mq. Alter- 
natively, we can estimate the Mhaio following the M^-Mhaio 
relation (Guo et a l. 2010), 

= n log \f M-^^o \ + ( ^-^^o V""" 

The fit is valid in the Mhaio range of lO^"-® to lO"-^ Mq, 
and thus describes better the most massive a. 40 galaxies. 
However, it cannot constrain the low mass Hl-bearing ha- 
los. As expected, such a M^-Mhaio relation, with a turnover 
Mhaio ^ 10^^ Mq, produces a similar trend of the (A) 
variation as seen in Figure [Mfb) for galaxies with A/* < 
2Q10.5 ]y[^ However, for the most massive galaxies with 
A/* > lO^"'^ Mq, such a trend disappears, and galaxies with 
similar M^, have similar A, regardless of their fni- Mean- 
while, (A) turns over to decrease with increasing Af* with 
approximately constant Mhi- 

As a result, the pattern of colors in Figure H^ b) predicts 
that the most gas rich ones with A/* near 10^° Mq should 
be associated with halos with the highest A parameters. We 
are testing this prediction by studying the stellar, gas and 
dark matter components of a sample of very high HI mass 
{Mhi > 10^" Mq), high gas fraction (for their stellar mass) 
galaxies, the HIghMass sample. Visual inspection of the 
HIghMass galaxies shows that the ones with intermediate M^ 
shows the strongest color gradient ( bluer in outer regions), as 
predicted for high A galaxies (He rnandez & Cervantes-Sodil 
i2006i) . Future work will constrain the spin parameters us- 
ing velocity fields derived from HI synthesis maps for 20 of 
the HIghMass galaxies. Furthermore, the HI distributions and 
sites of star formation will be examined to yield more detail 
on the SFL and possible mechanisms which inhibit the forma- 
tion of stars in these massive HI disks. 

The HI, stellar and SF properties of the ALFALFA pop- 
ulation as exemplified by the q;.40-SDSS-GALEX galaxies 
suggest that the Hl-selected population is biased towards ex- 
tended disks which reside in high A dark matter halos. As 
such, their disks are currently forming stars but in an ineffi- 
cient manner. Despite the low SFE, such galaxies can have 
comparable or even slightly higher current SFRs rel ative to 
the ones in the low A halos of similar halo mass (see ^ ^4.2.41 ). 
This combination of on-going star formation and inefficiency 
in the conversion of gas into stars causes their disks to be ex- 
tended, blue and of low metallicity, suggestive of their be- 
ing less evolved. Remembering too th at the Q!.40 galaxie s are 
among the least clustered population dMartin et alj|2012h . we 
are reminded that these systems likely follow a quite different, 
quiescent evolutionary history relative to galaxies residing in 
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higher density volumes. Although the details of the involved 
processes still elude us, the study of the ALFALFA popula- 
tion strongly suggests that where a galaxy is born and resides 
determines to a large extent how and when it converts its gas 
into stars. 

6. CONCLUSION 

Given the importance of proposed future extragalactic HI 
surveys at high redshift by the SKA, it is critical to develop a 
full understanding of the characteristics of gas rich galaxies at 
the present epoch. ALFALFA is an on-going blind HI survey, 
for the first time providing a full census of Hl-bearing objects 
over a cosmologically significant v olume of the local uni- 
verse. Building on the a. 40 catalog (iHaynes et al.ll20lTl) . we 
use the corresponding photometry available from the SDSS 
DR7 and GALEX GR6 catalogs to explore the population of 
galaxies detected by ALFALFA. The combined a.40-SDSS- 
GALEX sample includes 9417 galaxies. SED-fitting to the 
seven UV and optical bands yields the stellar properties of the 
Hl-selected galaxies, including their stellar masses, SFRs and 
internal extinctions. The lack of a correction for internal ex- 
tinction can lead to systematic underestimates of the optical 
luminosities and SFRs. In order to reduce the overestimate 
of internal extinction and thus the SFR with decreasing A/* 
while still accounting for the effect of dust in these mostly 
disk-dominated systems, we apply a prior distribution of ry 
in the fitting process. Although extinction is even more of 
an issue at the short wavelengths, the addition of photome- 
try in the UV bands is critical to the diagnostic power of the 
SFH because the UV measures break the degeneracy evident 
in optical-only colors, particularly among the red sequence 
galaxies. 

ALFALFA offers a complete statistical sampling of the full 
range of Mhi and Af,, from the most massive giant spirals 
with Mhi > 10^" M0 to the lowest HI mass dwarfs with 
log Mhi < 10^'^ M0, thereby providing a firmer basis for 
the derivation of the fundamental scaling relations linking the 
global properties. First, we confirm the existence of the re- 
lations which have been found in typical local optical-UV 
catalogs (e.g. ISalim et al.ll2007h . For example, (1) SFRs in- 
crease but SSFRs decrease with increasing M^,, with differ- 
ent slopes in the high and low M, ranges, with the transition 
occun-ing at A/, - 10^-^ M©; (2) the SSFR is tightly cor- 
related with the {NUV — r) color, especially after the lat- 
ter is corrected for internal extinction. Second, we investi- 
gate similar relations involving the HI mass. For example, (3) 
SFRs also increase but SSFRs decrease with increasing Mhi, 
though with larger scatter The HI gas contributes a signifi- 
cant fraction of the baryons in Hl-selected galaxies, and their 
SFRs show a good correlation with Mhi, suggesting that a 
global, atomic, volumetric SFL applies in Hl-selected sys- 
tems. (4) The HI fraction, fni = Mhi/M^,, nicely corre- 
lates with Af*, but a change in the slopes of the relation is 
evident near Af* ~ 10^ Mq. (5) Galaxies with bluer col- 
ors in general have higher fni- (6) The star formation ef- 
ficiency, SFE = SFR/ Mhi, mildly increases with stellar 
mass with a slightly steeper relation for M^, < 10^'^ M©. In 
we give the best linear fits to the principal scaling rela- 
tions among the ALFALFA population, including the Mhi- 
A/, and S S F R-M^ correlations, as well as the fundamental 
planes of fni-jNUV — r)-/x», etc. In particular, we argue 
that Equation (|2]|7]i provides the most robust predictor based 
on optical properties of the detection rate by future HI line 
surveys with the SKA and its pathfinders. 



Besides the scaling relations themselves, the combined 
a.40-SDSS-GALEX dataset, as a large and homogeneous 
sample with HI measures, enables the quantitative appraisal 
of the scatter in relations and a deeper understanding of the 
role HI plays in the galaxy evolution. In particular, the de- 
creasing fni with increasing A/» is related to the star-forming 
sequence identified in the SSFR-M^ diagram, or the evo- 
lutionary tracks of galaxies on the CMD, i.e. as galaxies 
assemble their stellar masses, they evolve gradually to rel- 
atively red and gas poor regimes, and also show lower SS- 
FRs. Furthermore, only evident among the low mass galaxies 
(A/, < 10^ ''^ Mq), the galaxies with higher fni on aver- 
age also have higher SSFR in a given A/, bin. Similarly, 
within a given Mr bin, higher fni on average indicates bluer 
(NUV — r) color However, the corresponding trend that the 
Hl-rich galaxies are more likely to be blue starburst galaxies 
with high SSFRs is weak among the high mass galaxies, i.e., 
the regulation of SF by HI is weaker in more massive sys- 
tems. Similarly, the dispersion of the color distribution in a 
given Mr bin and the dispersion of the SSFR distribution in 
a given Af* bin are both at a minimum near A/, ^ 10^'^ M0, 
and both increase at masses lower than that. It appears that 
SF is no longer regulated principally by the stellar mass in 
low mass systems. In their shallow potential wells, gas can be 
easily removed so that the SF quenches causing the galaxy to 
migrate onto the red sequence. 

We also focus on the nature of the population detected by 
the ALFALFA survey, in the context of populations better un- 
derstood through observations at other wavelengths, e.g. opti- 
cal or UV. The a. 40 galaxies on average have higher HI frac- 
tions in a given stellar mass bin, compared to the optically- 
selected samples, with an overall average of fni ^ 1-5. 
95.6% of the a.40-SDSS-GALEX galaxies have {NUV - 
r) < 4 and belong to the blue cloud on a UV-to-optical CMD. 
The red ALFALFA detections include early type galaxies with 
quenched SF but unusually high HI masses, suggesting a re- 
cent acquisition of HI. The very blue Hl-rich galaxies may be 
attributed to a SFH that steadily rises to the present day with 
little integrated past SF. The SFEs of the Hl-selected galax- 
ies are lower on average, or equivalently, their gas depletion 
timescales are longer (average tj^ ~ 8.9 Gyr), compared to 
the high stellar mass galax ies included in the GASS survey 
jSchiminovich et a"ni2010l) . Given the fact that the overall 
SFRs of a. 40 galaxies are even higher than those in GASS 
at a given Af,, the low SFEs found for Hl-selected galax- 
ies are caused by their high HI masses rather than by lower 
SFRs. This result is consistent with the idea that Hl-selected 
sample is biased towards the most gas-rich galaxies and that 
the SFE is low in Hl-dominated systems. A bottleneck may 
exist in the HI to H2 conversion, or the process of SF from 
H2 may obey an unusual SFL with low efficiency in the very 
Hl-dominated galaxies. For a given (FUV — NUV) color, 
Hl-selected galaxies have on average lower extinctions, sug- 
gesting that they have different SFHs or dust geometries. 

To quantify better the bias of the Hl-selected population rel- 
ative to the optically-selected galaxies, we constructed two 
volume-corrected control samples, starting from the a. 40 and 
the SDSS DR7 catalogs, which we designate Shi and Sopt re- 
spectively. The Hl-selected Shi sample is found to be biased 
against red-sequence galaxies as well as massive but low SFR, 
low SSFR galaxies. However, if only the blue cloud galaxies 
with {NUV — r) < 4 are considered, both samples define 
similar SF sequences, i.e., an HI survey well samples the star 
forming population. For the SDSS-selected volume-limited 
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Sopt sample, the rate of detection by ALFALFA decreases to- 
wards redder colors. Virtually all very blue Sopt galaxies at 
both the high and low stellar mass ends are detected by AL- 
FALFA; however, among the blue population, the HI detec- 
tion rate drops to ^ 40% throughout the high number den- 
sity region along the SF sequence. At the same time, only a 
very few of the galaxies which lie below the SF sequence in 
an SSFR vs. A/* diagram are detected by ALFALFA. Fur- 
thermore, ALFALFA misses a substantial fraction of the op- 
tical galaxies lying on the redder side of the blue cloud. The 
volume-corrected optically-selected Sopt sample well repro- 
duces various scalin g relations derived from the high stellar 
mass GASS sample ( Catinella et alj20TotlSchiminovich et alj 
I2OIOI) . However, at the highest stellar masses, the Hl-selected 
Shi galaxies show systematically larger discrepancies in their 
SF properties from the GASS results as the fraction of the to- 
tal population which is detected by ALFALFA in a given stel- 
lar mass bin declines. In comparison with optically-selected 
samples, Hl-selected galaxies that have high gas fractions 
are relatively less evolved and have, on average, bluer col- 
ors, higher SFRs and SSFRs, but lower SFEs, extinctions and 
metallicities. 

Previous authors, notably lBoissier & PrantzosI ( 120001) . have 
proposed that the overall low SFEs found in gas-rich systems 
may be explained if, their disks are characterized by large disk 
scale lengths and lower stellar surface densities because their 
host dark matter haloes have high spin parameters A. We ex- 
plore this hypothesis comparing the spin parameter distribu- 
tions of the volume-limit ed Shi and Sovt samples following 
the approach outlined in iHernandez et al.l (12007) which esti- 
ma tes the K-ot using the T ully-Fisher relation. As presented 
bv [Hernandez et alT(l2007 l). this estimate of A assumes that 
all galaxies are characterized by the same disk mass fraction 
Fd = 0.04. Under that (unlikely) assumption, we find a spin 
parameter distribution close to that found by Hernandez et al] 
^007) for their SDSS disk subsample and well-fit by a log 
normal distribution in agreement with numerical simulations. 
There is a slight hint that the Hl-selected population Shi has 
a slightly higher Ao than the Sopt sample which could re- 
flect the bias that the Hl-selected sample is characterized by 
more extended disks. However, abundance matching between 
the observed stellar mass functions and the CDM halo mass 
functions derived from simulations strongly suggests that that 
F is not a constant but r ather peaks arou nd a halo mass of 
Mhaio ~ 10^^ M© (e.g. lGuo et"ani20IOl) . At the low mass 
end, baryon depletion again grows in the shallow potential 
wells (Hoeftet al. 2006). Because the a. 40 catalog contains 
HI line widths, we calculate spin parameters from them adopt- 
ing the Vhaio-Vrot relation from Table 1 of Papastergis et al] 
(I2OI lb . While this method of assigning Mhaio is only vaUd 
in a statistical sense, the result is clear: the distribution of A 
is no longer log normal and has a mean value well in excess 
of the expectation, strongly reinforcing the hypothesis that the 
ALFALFA population favors high A dark matter hosts. 

In the future, our multiwavelength program to study the 
HIghMass sample of high gas fraction, high HI mass galax- 
ies will explore the star formation process and the interrela- 
tionships of the stellar, gas, dust and dark matter components 
within this set of exceptionally massive HI disks. The Survey 
of HI in Extremely Low mass Dwarfs (SHIELD) is explor- 
ing similar relationships among the lowest H I mass galaxies 
detected by ALFALFA dCannon etalj|2011h and has already 
uncovered evidence that the SFL in those galaxies is unusual. 
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Fig. 1 . — UV-optical color-color diagram. Contours and points represent the a.40-SDSS-GALEX common sample. The number density of galaxies in each 
gird cell is labeled on selected contours, e.g., the lowest contour level beyond which scatter points are plotted is 20 galaxies per grid cell. The grid size is shown 
by the interval of minor ticks on both axes. Among these 9417 galaxies, 96% have {NUV — r) < 4 and thus lie in the blue cloud by the criteria of ,Salim et alj 
<20Q7b . The {NUV - r) and (u - r) colors are well correlated among the blue galaxies, with a slope of 5{u — r)/S{NUV — r) ~ 0.6, but the distribution 
flattens among the reddest population. The (NUV — r) color serves as a stronger diagnostic of SFH than colors derived only from the optical bands. 
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Fig. 2. — Scaling relations between HI mass, stellar mass and color Contours and points follow the definitions in Figure[T] Blue diamonds and soUd lines 
indicate the (log y) values in each log x bin. The number of galaxies in each log x bin is Usted at the bottom of panels (c, d). Cyan dash-dotted lines in panels (c, 
d) denote the average values of the GASS galaxies with A/, > 10^" iCatinella et al.ll20T 9). Dashed lines in the same panels are from Cortese et al.l 1201 II) 
derived for galaxies belong to different environments: 'Hl-normal' (green on top), 'outside Virgo' (yellow in middle) and 'inside Virgo' (red at bottom). Typical 
error bars of individual galaxies are given in the comer of panels (a, d). Spearman's rank correlation coefficients, rs, are listed in all panels. 
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Fig. 3. — HI fraction estimators. Daslied lines are tlie one-to-one lines. Contours and points follow the definitions in Figure[T] Spearman's rank correlation 
coefficients are shown in the lower right comers of all panels. Upper panels - The ALFALFA-observed fui versus the predicted values based on the fundamental 
planes calibrated from the GASS sample (Catinella et al. 2010) in panel (a), as well as from the SDSS-selected sample (Zhang et al. 2009) in panels (b) and (c). 
Middle panels - Observed fuj versus the values predicted by the best fit to similar planes in this work, given the O.40-SDSS-GALEX sample. The systematic 
offsets are removed and the correlations are tighter. Bottom panels - Compared to the predictors in the middle row, the colors have been corrected for internal 
extinction and the stellar mass surface density is based on rgo.r, so that the scatter is reduced. 
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Fig. 4. — SF properties of the a.40-SDSS-GALEX population. Contours and points follow the definitions in Figure[T] Spearman's rank congelation coefficients 
are shown in the lower left corners of all panels. Colors here are not corrected for internal extinction. Typical en'or bars of individual galaxy estimates are plotted 
in the lower right corners of panels (a, e, i). The red dash-dotted line in panel (b), with a slope of 1.19, represents a global, atomic, volumetric SFL defined by 
the a.40-SDSS-GALEX galaxies. The red dash-dotted line in panel (d) represents the linear fit to the star-forming sequence given in equation (8). In the bottom 
row, tracing the SFE, the cyan dashed line shows the average value obtained bV iSchiminovich et al.i t2010i^ for the GASS sample, while the green dash-dotted 
line corresponds to the Hubble timescale. 
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Fig. 5. — Diagrams of the a.40-SDSS-GALEX sample with shade scale showing the averaged HI fraction in each giid, while number density is indicated by 
the contours. Left panel - UV-to-optical CMD. The blue cloud galaxies dominate in number and are associated with higher fni on average. At given Mr, a 
redder (NUV — r) color on average indicates lower fni', this trend is more evident at the faint end. Right panel - The star forming sequence as traced by the 
contours. As galaxies assemble Af* and evolve along the sequence, their HI fractions follow a decreasing trend. At a fixed M, , galaxies with lower SSFRs on 
average have lower /j//, which is also more evident among low mass galaxies with Af« < 10^'^ The broadenings of the distributions of color, SSFR and 
/j/j at the low mass end are correlated to each other. 
(A color version of this figure is available in the online journal.) 



a. 40: the Galaxy Population 



27 




Fig. 6. — Comparison of the basic properties relevant to the sample selection of the Hl-selected sample Shi ( 7157 galaxies selected from the a. 40, blue 
points and blue histograms) and the optically-selected one Sopt (16817 galaxies selected from the SDSS, red points and red histograms). The panels above the 
histograms show the fraction of the Sopt galaxies that are cross-matched to a. 40 in each bin, similar to an HI detection rate of the Sopt galaxies by ALFALFA. 
Both samples are extracted from the same sky area (panel c), lie within cz = 15000 Itm (panel d) and are cross-matched to the GALEX catalog. Spaenhauer 
diagrams are in panels (a) and (b), i.e. r-band absolute magnitude, Mr, and M[j;i against distance. A weight (or volume correction) cut of 60 applied on Sopt 
results in Mr brighter than ~ —16, whereas applied on Shj results in Algj > 10*'^ Mq. Though the two samples probe similar Mr ranges (panel e), the 
SDSS is deeper than ALFALFA, as evident in the distance distribution (panels d). No bias against edge-on galaxies in the Shi population is seen in panel (f). 
(A color version of this figure is available in the online journal.) 
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Fig. 7. — Panels (a-c) - The r-band internal extinction versus stellar mass. The values in panel (a) are derived by SED-fitting without the application of a 
prior ty distribution. The mean of the prior distribution is plotted in panel (b), based on equation (12) in 'Giovanelli et al] 1 19971) . The SED-fitting has been 
improved by adopting the prior Ty distribution as evident in panel (c). Internal extinction is a weakly increasing function of M* with the greatest correlation, 
with a coefficient = 0.33, shown in panel (c). Panel (d) - Internal extinction is systematically higher in more inclined galaxies, as expected. 
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Fig. 8. — SSFR versus intrinsic colors, denoted by tlie subscript zero, after correction for internal extinction. The left column shows results for Sopt while the 
light shows Shi- Pearson congelation coefficients and typical error bars are also shown and the median intrinsic colors are indicated. In addition to the blueward 
shift of the distribution in panel (b) the scatter here is also greatly reduced compared to the similar plot shown in FigureOf) where no internal extinction correction 
has been applied. The (NUV — r) color breaks down the degeneracy of (u — r) in the red range when inferring the SFH. The red and low SSFR tail in Sopt 
disappears in Shi- Shi is on average intrinsically bluer than Sopt in both colors. 
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Fig. 9. — Computed Apuy - UV color diagrams for Sopt in the left column and Shi in the right column. The upper row shows the results for the galaxies 
themselves, without weighting while the bottom row shows results after applying the Vsur /Vmax weighting scheme. Typical error bars are shown in the corners 
of panels (a) and (b), together with the median Apjjv values listed for both the samples. The red dashed line coiresponds to the fit to the star-forming galaxies 
derived by Salim et al. (2007), based on a typical local SDSS-GALEX cross-matched catalog. It is in agreement with Sopt- In contrast, the Shi distribution 
is offset from the fit, with lower Apuy at a fixed {FUV — NUV) color, which may due to lower metallicity, different SFH, and/or dust geometry of the 
Hl-selected galaxies. 
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Fig. 10. — Intrinsic optical color-magnitude diagra m of Sppt (left pane ls) and Shi (right panels); before (upper) and after (lower) applying the Vsur /Vmax 
weight correction. The dash-dotted curve is derived in' Baldrv et alj (20040 as the best fit to the division of red sequence and blue cloud, shifted for the extinction 
connections. Typical error bars including the uncertainties in the extinction corrections are plotted in the lower left comer of panels (a) and (b). The bimodel 
distribution is more evident in the Sopt representations; colors become generally bluer in fainter galaxies for both populations. The weighted panels better 
represent the luminosity function, which predicts more faint galaxies relative to the bright ones in the blue cloud. However, a second peak of number density 
at the faint end is seen on the red sequence, which suggests a diversity of SFH within the dwarf galaxies. A higher percentage of galaxies lie bluewards of the 
division in Shi (84%) relative to Sopt (67.5%). HI surveys such as ALFALFA are highly biased against red sequence galaxies. 
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Fig. 1 1 . — Diagrams with sliade scale sliowing the fraction of galaxies in Sopt that are cross-matched to the a. 40 in each grid, which is close to the HI detection 
rate of the Sopt by a AO. The Sopt number density is indicated by the contours. The cross-match fraction is the lower limit of the detection rate of Sopt galaxies 
by 0.40; its overall average is 34%. Left panel - Intrinsic UV-to-optical CMD of Sopt- ALFALFA is more efficient in detecting blue galaxies, especially (i) the 
very bright and blue galaxies with huge gas reservoirs and (ii) the galaxies with the highest HI fraction lie on the faint end of the blue cloud. However, starting 
from the redder edge of the blue cloud, the cross-match rate drops below 10% and to even ~0% throughout most of the red sequence. Right panel - SSFR versus 
stellar mass for the Sopt galaxies. The cross-match rate is the highest among the galaxies with high SSFRs at both high and low stellar mass ends and is close 
to the overall average thi'oughout the high number density region along the star forming sequence. However, it drops to below ^^20% from the lower edge of the 
sequence to even ~0% in the low SSFR regime. 
(A color version of this figure is available in the online journal.) 
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Fig. 12. — SFR versus M« for the Sopt (left) and the Shi (right) samples; upper panels show the individual galaxies while the low ones show the results of 
applying the Vaur /Vmax weighting. The cyan dashed Unes above A/» = 10^^ Mq show the fit to this relation obtained by Schiminovich et alj <201(]|) . based 
on the GASS sample. The blue diamonds and lines represent the corresponding log(S'-F-R) of our datasets. The number of galaxies in each stellar mass bin is 
listed at the top in panels (a) and (b), together with the typical error bars plotted in the lower light corners. The median SFRs for both the samples are indicated. 
The GASS fit is consistent with the distribution seen for the Sopt sample but is systematically below the average of the Hl-selected galaxies. The slope of the 
relation appears to steepen below A/* 10^'^ Af©. The Shi sample has a higher overall SFR value than the Sopt galaxies, and probes to sHghtly lower A/* 
ranges with generally lower SFRs. The data points reveal a concentration of massive low SFR galaxies below the main trend in Sopt, which are largely absent 
from the Shi population. An HI survey samples the star-forming population. The flattening of the trend in the lowest mass bin in the Sopt plots is artificial, due 
to the weight cut applied. 
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Fig . 13. — Similar diagrams of SSFR versus M« . The numbers of galaxies in each stellar mass bin are listed, as well as the median SSFRs for both samples. 
The contours outlying the high n umber dens ity region roughly trace the star forming sequence. The red dashed hne shows the fit to such a sequence of the blue 
galaxies with {NUV — r) < 4 (Salim et"a l. 2007), derived from a typical local SDSS-GALEX cross-matched catalog, with the majority confined to the stellar 
mass range 10^ — IQ-'^'^ M (7). The blue diamonds and lines are obtained by applying the same color criteria to the Sopt or Shi galaxies. The cyan dashed line 
comes from lSchiminovicF et al. (2010) based on the high A/, GASS sample; it well represents the contours of Sopt, but lies systematically below the average 
of Shi- When only blue galaxies are considered in both samples, the discrepancy between the main trend of Sopt and Shi is small. Furthermore, both agree 
well with the red dashed line. Galaxies selected by HI criteria have on average higher SSFRs than optically-selected ones. The breakdown of the star forming 
sequence above A/, ~ 2 X 10^^ Mq is only evident among the Sopt sample, whereas the Shi galaxies are strongly biased against massive and low SSFR 
galaxies. 
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Fig. 14. — Spin parameter A distribution obtained by assigning tlie M^alo tlirough a Vrot-Vhalo relationsliip. Panel (a) - The black solid line is the normalized 
PDF of the A distribution of 7459 galaxies from the O.40-SDSS-GALEX sample. The blue dashed line is the best lognormal fit to the distribution. The red dash- 
dotted line is the lognormal fit to a sample of spiral galaxies from the SDSS, assuming constant baryonic mass fraction F = 0.04 (Aq = 0.0585 Hernandez et aj] 
120071) . The A distribution of the a. 40 galaxies has a higher mean A (0.0852) and a wider dispersion, mainly arising from the adopted distribution of baryonic 
fraction. Panel (b) - fui vs. Af» diagram color coded by the mean A values of the galaxies in each grid. Below A/» ^ lO^'''^ M©, (A) increases with increasing 
fuj at a given A/*, but remains almost constant along constant Mhi lines. Such a trend reflects the adopted variation of F: Hl-selected gas-rich galaxies favor 
high A halos. 

(A color version of this figure is available in the online journal.) 



